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Abstract. We obtain Plancherel-Rotach type asymptotics valid in all regions of the complex plane for 
orthogonal polynomials with varying weights of the form e ~ N V( x ) on the real line, assuming that V has only 
two Lipschitz continuous derivatives and that the corresponding equilibrium measure has typical support 
properties. As an application we extend the universality class for bulk and edge asymptotics of eigenvalue 
statistics in unitary invariant Hermitian random matrix theory. Our methodology involves developing a 
new technique of asymptotic analysis for matrix Riemann-Hilbert problems with nonanalytic jump matrices 
suitable for analyzing such problems even near transition points where the solution changes from oscillatory 
to exponential behavior. 



1. Introduction 

Let V(x) be a real-valued function (an external field or potential) growing faster as \x\ — > oo than 
[log(l + x 2 )] 1+e for some e > 0. In particular, this implies that all moments of the measure on K given by 

(1) dv N {x) := e - Ny W dx 

are finite. A measure of this form is said to be a varying weight due to the presence of the parameter N . This 
paper concerns the asymptotic behavior of polynomials orthogonal with respect to varying weights of the 
form ([!]). They are defined as follows. For n = 0, 1,2, . . . , let p n — p n {z) = p n (z; N) = n n , n z n + • • • + k„ i0 , 
K n n > be the (unique) polynomial of degree n satisfying 



(2) / p n {x)x K dv N {x) = , for < k < 

Jm 

and 

(3) / p n (x) 2 dv N (x) = 1. 



The interest is in the behavior of the polynomials of degree N and N — 1 , where the integer N is the same 
which appears in the measure of orthogonality, in the limit N — > oo. We obtain a precise description of 
the polynomials pn(z; N) and pn-i(z; N) which is uniformly valid for all z £ C, for all N sufficiently large. 
This type of asymptotic description is often referred to as Plancherel-Rotach asymptotics, after the analysis 
of the Hermite polynomials in [21] . 

In the late 1990s new Riemann-Hilbert techniques originally developed for the asymptotic analysis of 
problems in integrablc nonlinear partial differential equations were applied to the asymptotic analysis of 
Riemann-Hilbert problems encoding systems of orthogonal polynomials with respect to varying weights on 
K, first for the case of external fields of the form V(x) — x 4 — jx 2 [2], and then for the case of general real 
analytic V in |6j[7j. (See [6] for more information about Plancherel-Rotach type asymptotics for orthogo- 
nal polynomials prior to the use of Riemann-Hilbert techniques.) The Riemann-Hilbert method has been 
extended, and applied to various types of asymptotic questions in approximation theory (see, for example, 
[12"] . [T5] , [T], and [T3]). With the exception of [T2] and [TS], all of these applications and extensions deal 
with weights that are real analytic. 

The main purpose of this manuscript is to establish Plancherel-Rotach type asymptotics for orthogonal 
polynomials, when the external field V possesses only two Lipschitz continuous derivatives, i.e. in the absence 
of analyticity. (The precise assumptions on the external field V are most naturally described in terms of the 



equilibrium measure to be defined in subsection 1.3 below.) To obtain such a uniform asymptotic description 
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we present a new hybrid Riemann-Hiibcrt-<9 method of asymptotic analysis, that is a significant extension of 
the d- method introduced in |15j to analyze orthogonal polynomials on the unit circle. By contrast with that 
method, a fundamental new feature of orthogonal polynomials with varying weights on the real line is the 
presence of "transition points" (also known as endpoints of the equilibrium measure) in the neighborhood of 
which the asymptotic behavior exhibits a complicated transition from oscillatory to exponential behavior. 

1.1. Application to random matrices. Among many applications of the asymptotic theory of orthogonal 
polynomials is the calculation of certain statistics of eigenvalues in random matrix theory. Unitary invariant 
ensembles of random matrices are described by probability measures of the form 

(4) dP N (M) = J- e -NTr(V(M)) rfMj 

Zn 

defined on N x N Hermitian matrices M, where V(x) is an external field of the type described earlier. Here 
e?M denotes Lebesgue measure on the algebraically independent entries: 

N 

(5) dM=Y[dM jj Yl dRe(M jk )dlm(M jk ), 

j=l l<j<k<N 

and Zn is a normalization constant (partition function) . One of the origins of the theory of random matrices 
in physics was the study of nuclear resonance levels in the 1950s. See [TB] and the references contained 
therein for more information. 

1.1.1. Connection with orthogonal polynomials. A remarkable connection to orthogonal polynomials was 
discovered in the late 1960s by Gaudin and Mehta [17] • The connection is the following formula for the 
density of the probability measure on eigenvalues induced by Q : 

(6) F N (\ 1 ,...,\ N ) = -^det(K N (\ i ,\ j )) 1 < ij < N , 

where the function K^{x 1 y) is the so-called reproducing kernel of orthogonal polynomials: 

N-l 

(7) K N (x,y) = e -^(v(x)+v(v))/2 £ Pnix)pn{y)> 

the polynomials p n {x) being defined in ([2])-([3]). It is a basic result of the theory that ^ indeed defines a 
probability measure on R . 

From formula Q one may effectively compute many statistical quantities involving the eigenvalues (see 
[IB] , and also [3]). Two examples are as follows: 

• Mean density of eigenvalues p\ N (A) defined as 

(8) p[ N \x) := ~ E N ( {eigenvalues Xj such that Xj < A} 

dX V JS 

where Ejv(-) denotes the expectation of • with respect to the probability measure ^ or equivalently 
This may be equivalently represented in terms of the orthogonal polynomials: 

(9) p[ N \x)=^K N (X,X)- 

• Gap probabilities F( a b ^ defined as 

(10) ^(a.h) : = Prob (no eigenvalues in (a, b)) , 

which may be equivalently represented in terms of a Fredholm determinant built out of the orthogonal 
polynomials: 

(11) P{a,b) = det(l-lC N \v (a , b) ). 

Here the integral operator JCn ■ L 2 (a,b) — * L 2 (a,b) possesses the integral kernel Kn(x,u): 

(12) JC N h{x)= [ K N (x,y)h(y)dy. 



One important example of the gap probability described in (10 1 and (11 1 is the case that b = oo, for then 
the gap probability coincides with the distribution function of the largest eigenvalue: 

(!3) F {a:+oo) = Prob (A max < a) = det (l - IC N \ LHa>+oo) J . 

1.1.2. Asymptotic behavior as N — > oo. A basic and important result concerning the N — > oo asymptotic 
behavior of random matrices is that the mean density of eigenvalues p[ N ^ has a limit: for all A G R, 

(14) Jim p[ N \x) = -0(A). 
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Note: the Gaussian Unitary Ensemble (GUE) first studied by Wigner corresponds to V(x) = x 2 , and in this 



case "0(A) = it x \j2 — A 2 , which is the famous Wigner semicircle law. It is well-known that the limit (14 1 



exists for quite general V(x). It is also known that if V is real analytic, the convergence in (14 1 is uniform. 
For those nonanalytic V for which existence of the limiting density ip(X) can be established, the convergence 
implied by the statement (|14|) has only been proven in a weaker sense. One consequence of the present work 



is that the convergence in ( 14 1 is in fact uniform assuming only that the function V possesses 2 Lipschitz 
continuous derivatives. 

The function tp is also a well-known quantity in approximation theory, where it is referred to as the density 



of the equilibrium measure. The equilibrium measure is defined in subsection 1.3 (for the purposes of the 
current discussion one may take the parameter c appearing in the definition of the equilibrium measure to 
be unity). 

In many circumstances, the largest eigenvalue distribution has been shown to possess a limit as N — > oo 
known as the Tracy- Widom distribution, a distribution function expressible in closed form in terms of the 
Hastings-McLeod solution of the Painleve II equation. The form of the asymptotic result is: 

(15) lim Prob (A max < (3 + (AA0~ 2/3 s) = F TW (s) 

where the constant A depends on the external field V, (3 — sup (supp( ?/>)), and i 7 Tw(s) is the famous Tracy- 
Widom distribution. 

Another fundamental object concerning the eigenvalues of random matrices is the limiting spacing distri- 
bution. Properly speaking, this is defined in terms of the spacing between ordered eigenvalues; however a 
"poor-man's" version of this distribution is the following (easier to define) quantity: 

(16) Q(s) := lim Prob no eigenvalues in la, a + 



N ^°° V V ' Np[ N \a) J J ' 

This limit is known to exist provided the external field is real analytic and provided that a is such that 
0(a) > 0, and it turns out that the function Q(s) which emerges in the limit is universal in that it does not 
depend on properties of the function V. Indeed, under the assumption that V is real analytic, one has 



(17) Q(s) = det [1 - S\^ Qt , : 

where S is an integral operator on the interval (0,s): 

„, , s f s sin(7r(.T — y)) , . . , 

(18) Sh(x):=J o ^_f h(y)dy. 



Via the connection to orthogonal polynomials explained earlier, the following asymptotic result concerning 
the reproducing kernel K^(x,y) built from the orthogonal polynomials implies (15 1: 

Asymptotic Result 1. There is a constant A so that for every u,v G 



(19) 



1 



N^oc (A N) 2 / 3 



K N [(3 + 



we have 
Ai(u)Ai'(«) - Ai(«)Ai'(u) 



(AA) 2 / 3 '^ (AAO 2 / 3 

Here Ai denotes the unique solution to Airy's equation y" = xy that is real, and that behaves as follows as 
x -» +oo: Ai(ar) - e~ 2xV2 / 3 /(2^x 1 / i ). 



Similarly, the limit appearing in (16 1 is implied by the following result: 
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Asymptotic Result 2. For every a with ip(a) > 0, and every u,v £ K, we have 

1 -rr ( u v \ sin(7r(w — v)) 

(20) lim K N [a+ a-' 1 - 



Asymptotic Result [l] was first established in the special case of the Gaussian Unitary Ensemble (i.e. 
V(x) = x 2 ) [22], using the classical Plancherel-Rotach asymptotics of Hermite polynomials [21]. This 
was extended to quartic potentials of the form V(x) — x A — 7a; 2 in [2], where furthermore Asymptotic 
Result [2] was also established. Because the polynomials associated with quartic V are not known to possess 
elementary contour integral representations, the analysis of [2 required a new method, namely the use of the 
Riemann-Hilbert formulation of orthogonal polynomials found in [TU]. Asymptotic Result [2] was established 
for general real analytic potentials V in [7], and the asymptotic formulae for orthogonal polynomials given 
in [7] were used to establish Asymptotic Result [l] in [5]. The new strategy introduced in [7] was a general 
method linking the equilibrium measure associated with V to a so-called g-function enabling the use of the 
non-commutative steepest descent technique for Riemann-Hilbert problems originally invented by Deift and 
Zhou [SJ and extended in [8 . Pastur and Shcherbina [T!J] have also studied the problem of establishing 
Asymptotic Result [2] under the assumption that V has three continuous derivatives. 

As is clear from the above discussion, the historical trend is toward establishing Asymptotic Results [l] 
and [2] for more and more general external fields V. The program of universality in random matrix theory 
is concerned with determining the most general external fields V under which such results hold true. In 
particular, it is of some interest to admit external fields that are not real analytic. As pointed out by Deift in 
[1] , the steepest descent method that works so well for analytic V cannot be easily applied to the nonanalytic 
case. The authors recently introduced a "9 steepest descent method" applicable to some Riemann-Hilbert 
problems involving nonanalytic data, but as formulated in |15j this method does not apply to the orthogonal 
polynomials described by conditions ([2])-([3]) because the equilibrium measure is compactly supported and 
the endpoints of support obstruct the type of nonanalytic deformations involved in the method. 

Among the applications of the results in this manuscript are proofs of Asymptotic Results [l] and [2] under 
weakened hypotheses on the external field V (we require two Lipschitz continuous derivatives) via rigorous 
Plancherel-Rotach type asymptotics for orthogonal polynomials. Our method involves a hybrid "Ricmann- 
Hilbert-c* steepest descent method" generalizing the simpler method of [TS] to handle support endpoints. 
Our results hold under weaker conditions on V than those under which Asymptotic Result [2] is considered in 
the paper [T^] , and to our knowledge we have the first proof that Asymptotic Result [T] holds in the absence 
of analyticity of V. 

We remark that it is not necessary to first obtain asymptotics for the orthogonal polynomials themselves 
in order to deduce enough information about the reproducing kernels Kj^(x,y) to establish Asymptotic 
Results [l] and [2] for certain general external fields V. For example, a recently introduced method (based on a 
comparison principle for Christoffel functions) of Levin and Lubinsky has been quite successful in establishing 
Asymptotic Result [2] [14] under extremely weak global conditions on the external field V and stronger local 
conditions (but still far from analyticity) near the point a of expansion in the spectrum. Also, Asymptotic 
Result [T] has been studied for real analytic V without the use of orthogonal polynomials by Pastur and 
Shcherbina [201- 



1.2. Essence of the d method. In the asymptotic analysis of Riemann-Hilbert problems there is an analog 
of contour deformations which plays a crucial role in identifying subsets of the plane which produce the 
dominant contribution to the Riemann-Hilbert problem's solution. It is common to begin with a Riemann- 
Hilbert problem whose solution, A, is analytic off a given contour £a, and to use explicit piecewise analytic 
quantities to define a new matrix B solving a new equivalent Riemann-Hilbert problem in which the relevant 
contour Sb is a deformation of the original contour £a- 

A fundamental obstacle to this procedure occurs when one requires the analytic extension from a given 
contour of a rapidly oscillating function whose phase possesses no analyticity properties. For the asymptotic 
analysis of orthogonal polynomials with varying weights on the real line, in which the external field V 
possesses only finitely many derivatives, this is a central issue. 
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In this paper, the new approach which circumvents this problem is to depart from Riemann-Hilbcrt 
problems entirely, by introducing transformations that explicitly violate analyticity. Instead of Riemann- 
Hilbert problems, we then characterize our newly-defined matrix- valued function as the unique solution of a 
d problem. 

Given a smooth matrix-valued function W(x, y) of compact support in K 2 , a d problem is a first-order 
system of linear partial differential equations on M. 2 involving W(x, y) as coefficients and the Cauchy-Riemann 
operator 

acting on the unknown. Here is a prototypical example. 

d Problem (Prototype). Determine a 2x2 matrix A(x, y) for (x, y) G M 2 having the following properties: 
Continuity. A(x, y) is a continuous function of x and y for x + iy G C. 
Deviation From Analyticity. For x + iy G C. 

(22) dA(x,y)=A(x,y)W(x,y), 

(note that A(x, y) is analytic off the support o/W, because there one has <9A(x, y) = 0). 
Normalization. The matrix A(x,y) is normalized as follows: 

(23) lim A(x,y)=I. 

x.y — >oo 

Once one admits the possibility to use non-analytic extensions of functions originally defined on contours, 
one is faced with an overabundance of choices, and the issue becomes one of selecting, constructing, or 
otherwise establishing the existence of, an extension suitable for subsequent asymptotic analysis. 

This idea actually yields an interesting approach to a classical result of asymptotic analysis. Given a 
real- valued function 9 : [— 1, 1] — > M satisfying 9(0) — 9'(0) — 0, 9"(0) > 0, the problem is to provide a large 
n asymptotic description for the integral 

(24) 7(n) := J e in9 ^ dx . 

For convenience, let us assume that 8"'(x) is bounded and 9"(x) > w > for all x G [—1, 1]. 

The usual approach to this problem (see, for example, [18) involves many steps, including integration by 
parts, implicit variable changes, and Taylor expansions, the result of which is 



(25) 'M = JJk> e '""( 1 + ° i."' 1 ' 2 ))- 

We may instead establish this in the following way, which elucidates certain aspects of the methods we use 
in the sequel. Let 0(x, y) represent an arbitrary extension of 9(x), which satisfies 0(x, 0) = 8(x). Then with 
the aid of Stokes' theorem, we may write 

(26) 7(n) = - J e me( - x ^ d(x + iy) + 2i J J de me( - x > y) dA, 

where T represents a contour in C from 1 to —1 (different than the interval [—1,1]), and A represents the 
(oriented) area enclosed by the oriented contour formed by the union of [—1, 1] with T. See Figu re [l] 

Were the function Q(x,y) analytic, the double integral appearing on the right-hand side of (|26[) would 
not be present, and we could choose the contour V to be the contour of steepest descent. Although <d(x,y) 
cannot be analytic if only three derivatives of 9(x) are assumed to exist, we nonetheless observe that the 



right-hand side of ( 26 ) is still independent of the choice of both the contour T and the particular extension 



0(x, y). This begs the question: Can we pick the extension 0(x, y) and the contour T so that the right-hand 



side of ( 26 1 may be easily estimated? The answer is yes. 

We take T to be the contour comprised of a vertical segment from (1,0) to (1, 1), followed by the line 
segment connecting (1, 1) to (—1,-1), and ending with the vertical line segment from (— 1, —1) to (—1,0), 
and let A + and A_ denote the interior of the two triangles formed by this contour and the real interval 
[—1,1], A + in the first quadrant, and A_ in the third quadrant, exactly as illustrated in Figure [T| We 
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will explicitly construct an extension 0(x,y) of 0(x) defined for (x,y) G A + U A_ to satisfy the following 
conditions for some constants K > and k > 0: 

(CI) 6(x,0) = 6(x), for -1 < x < 1. 

(C2) 9(x,x) = |6<"(0)(a; + ix) 2 , for -1 < x < 1. 

(C3) |<96(x, y) I < i^Tj/ 2 for all (x, y) € A+ U A_. 



(C4) Im(6(x, y)) > kxy for all (x, y) G ^4+ U A_ 



Using such an extension and properties (CI) and (C2), the representation (26) may be rewritten as 



\/2 

(27) I(n) - e"/ 4 / e ~ nf> "^ s ^ 2 ds 

-V2 



i I e 
'o 



1 ^&{-x,y) dy + l f° e *«o(i,») 



dy 



2n / / e in@ ^dQ{x, y) dxdy + 2n / / e^^defo y) dx dy . 

JJa + J J A- 



The four integrals on the right-hand side may be estimated directly with the help of properties (C3) and 
(C4): 



(28) 



% I e 

o 



mO( — 



dy 



I e - n Im(&(-l,- S )) ds < /" e -to ds < 

o Jo 



e- fens ds 



1 

fcn 



(29) 



i / ' 



m0(l,y) 



: / e -«M©(M) ds < / e -fc«s ds < , 
'o Jo Jo 
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-*™ ds = — 
kn 



%e ^de(x,y)dxdy 



< 2n 



-nlm(0(a:,y)) 



< 



(30) 



dQ(x, y)\ dxdy 
2Kn J J e- knxy y 2 dxdy 
2Kn J J e' knxy y 2 dxdy 



< 2Kn 

_ 2K 
k 2 n 



■ dr 



7r/4 



-U „,a 



^ e -fe„r 2 co S (e) sill (e) r 2 sill 2^) 

/ 4 de 

cos 2 (<9) ' 



Therefore all terms on the right-hand side of (27) are 0(n *) as n — > oo. Now, since 

-V2 



(31) 



t/4 



e -«e"(o) s - ! /2 rfs = 



2tt 



ra0"(O) 



r/4 



exponentially small terms as n — > oo, 



we have established ( |25| ) if we can find an extension 0(a;,?/) of 9{x) satisfying conditions (C1)-(C4). 

The extension <d(x,y) may be defined as follows. First let B(t) represent a "cut-off" or "bump" function 
which is infinitely diffcrentiable and satisfies B{t) = for t near 0, and B(t) = 1 for t near 1. More precisely, 
we assume that B : R — ► [0, 1] is of class C^°°^(IR) and satisfies B(x) = for x < and B(x) = 1 for x > 1. 
An example of such a function is 



(32) 



B(t) := { - tanh 



t 



1 -< 2 



1 

2 ' 



x < 
< x < 1 



1, x > 1, 

but our analysis will never require the detail of this formula. Next define 



(33) 



e (x,y) := 6{x) + iy6'(x) + l -(iy) 2 6" (x) , 
e*°\x,y) := l -6"{0){x + iy) 2 . 



Our extension Q(x,y) is then defined via 



(34) 



0(x,y) 



1 — B 



Q (x,y) + B(^)e^\x, y ). 



Note that the function <d (x,y) is an extension of 9(x) that satisfies 



(35) 



de (x,y) = ^(iy) 2 6'"(x). 



The function Qo(x,y) is a rectilinear version of the type of extension discussed in [15j . It does not match 
the desired quadratic on the diagonal part of the contour T; we use the function B to smoothly deform this 
extension to the quadratic Qq o1 (x, y). Straightforward calculations show that 0(x, y) defined in (34 1 satisifes 
the four conditions (C1)-(C4) described above. Indeed, 6(x,0) = 0o(x, 0) = 9{x), so condition (CI) holds, 
and Q(x,x) = Oo lol (a;,a;) = h0"(0)(x + ix) 2 , so condition (C2) holds. To confirm condition (C3), note first 
that 3Qq° 1 (x, y) = 0, so using (35 1 we have 



1 - B 



(!) 



(36) dQ(x,y) = - 
and then by Taylor expansion 

(37) e (x,y)-e*°\x,y) 



le"<{x)y 2 + -B< 
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(I) 


' y .i" 




X z X 




"\t2)x 2 y 



[e (x,y)-G*°\x,y)] , 



'(^)xy 2 , 



for some numbers £1, £2, and £3 in [—1, 1], so since 6"' is uniformly bounded and < y/x < 1 for (x,y) € 

eo(x,y)-e* ol (x,y)=0(x 3 ), 



A + U , we have 
(38) 
and so 
(39) 



d&{x,y)=0{y 2 



B'(^jO(x 2 ). 



Finally, since B'{t) = 0(t 2 ) holds for all t £ R, we have confirmed condition (C3). To check condition (C4), 
we calculate directly 



(40) 



Im(6(x, y)) = [l - B (I)] lm(G {x, y)) + B (|) Im(6|| ol (x, y)) 

y6'{x) + B (|) 0"(O>y 
^"(0 + s(|)^(0)a:y, 



1 -B 
1 -5 



for some number £ £ [—1,1]. Then since by assumption 8"(x) > w > holds for x € [— 1,1] and since 
B [0,1], we have 



(41) 



Im(6(x,y)) > 



1 - B 



(!) 



uixy = wxy, 



so condition (C4) is verified as well. 

Now we will be starting with a 2 x 2 matrix B, which is the solution of a Riemann-Hilbert problem in 
which the jump matrix contains entries of the form e ln9 ( x ) ^ with 9 real, and possessing only 2 Lipschitz 
continuous derivatives. We will define an extension of 9 in exactly the spirit of the above calculations, and 
use it to define a new matrix-valued function D, that is no longer analytic. The matrix-valued function D 
will be characterized by a hybrid Riemann-Hilbert- d problem. The main point is this: our extension of 9 
will be chosen so that this hybrid Ricmann-Hilbcrt-9 problem succumbs to a large-n asymptotic analysis. 

1.3. The equilibrium measure and associated quantities. The so-called equilibrium measure associ- 
ated with the function V(x) and the ratio c := N/n is well-known to be a key ingredient in large-degree 
asymptotics of the polynomial p n {z) of degree n in the orthonormal system associated with the measure 
duisr(x) defined in terms of N and V by 0. Here c > is held fixed as n (and hence also N) tends to 
infinity. Generally, given a real-valued field V(x) defined for x £ R and a parameter c > 0, we may consider 
the following associated weighted energy of a positive charge (Borel measure) [i on the real line M C C: 



(42) 



log 



PP(m) "'suppO) 



>l\ 



dji(x) dfx(y) 



V(x) d^{x). 



suppO) 



The equilibrium measure /Lt* is defined to be the unique positive measure ^ minimizing E[fi] subject to the 
constraint 



(43) 



d/j,(x) = 1. 



'suppO) 

The equilibrium measure is equivalently characterized by the corresponding Euler-Lagrange variational con- 
ditions. There is a real constant £ (the Lagrange multiplier originating from the constraint (43)) such that 

1 



(44) 

and 
(45) 



SE 
6fi 



SE 
5(i 



(x) := 2 



(x) 



log 



supp(^») 



log 



suppO.) 



\x - y\ 



F - y\ 



d^*{y) + cV(x) = — £, x £ supp(/x»), 



d^*(y) + cV(x) > — £, x $ supp(/z*). 



1.4. Assumptions on external field V. We now impose several conditions on the external field, some of 
which are best described in terms of the equilibrium measure and its complex valued "log-transform" g(z) 
defined below in ( 46 1 . 

Condition (Smoothness of V). The external field V possesses two Lipschitz continuous derivatives. 

A consequence of this is that the equilibrium measure is absolutely continuous with respect to Lebesgue 
measure, with continuous density if>{x). 

Condition 1 (Support properties of //*). We suppose that the external field V is such that the equilibrium 
measure is supported on a finite union of intevals, U^-J^o,-, /3j], with ax < fix < 012 < P2 < • • • < Pg+i- 

By convention for future convenience, we set /3o := —00 and ac+2 '■= +00. 

To describe further conditions on V imposed via its equilibrium measure //* , we will require an auxiliary 
function g(z) analytic for z £ C \ (— oo,/?g+i) defined in terms of /j* by 



(46) 



g(z) := / log(z -s)d^(s) 



log(z — s)ip(s) ds, 



where ip(x) is the Radon-Nikodym derivative of fj,*, that is, dfj,*(x) = tp(x)dx. Here we are choosing the 
branch cut of the integrand so that for each s£l, log (z — s) is an analytic function of z for z £ C \ (—00, s] 
that is real- valued for z > s, which ensures the claimed analyticity properties of g{z). In terms of g(z) the 
variational condition ( |44| becomes 

(47) cV{x) - (g+{x) + g-(x)) = -i, ay < x < fa, j = l,...,G+l, 

where g+(x) and g~{x) denote the boundary values taken by g(z) as z — ► x with z € C+ and z G (Ir- 
respectively. Also, Condition 1 and the reality of the equilibrium measure together imply that there are real 
constants £Iq, . . . , flc such that 



(48) 



iQ , 

9+( x ) ~9-{ x ) = { 

0, 



x < ax, 

(3j < x < ay+i , 
x > (3g+i , 



J 



1, 



,G, 



and from the normalization (43 1 it follows further that f2o — 27r. Furthermore, since /i* is a positive measure, 
(49) 9{x) :=-i( 9+ (x) - , g_[x)) 

is real and nonincreasing for all x € R, so in particular 2ir = il > fij > • • • > > 0. Assuming that 
differentiation commutes with taking boundary values (this may be easily justified later) ( 47 1 and ( 48 1 imply 
that 



(50) 



9+{x) + g'-(x) = cV'(x) , aj < x < fy, j = 1, 
g' (x) - g'_ (x) = 0, x € K \ supp(^) . 



,G+1, 



In particular, g'(z) is an analytic function for z £ C \ supp(V')- 
Next, for x £ K, define the real-valued function 



(51) 



(f>(x) := cV{x) +l-g+(x)- g-(x). 



According to (44 1 and (45), we have <f>(x) = for x £ supp(?/>) and <fi(x) > for \ supp(^). 

Finally, for j = 1, . . . , G + 1 we define functions h aj : - 
formulae 



(52) 



, G + 1 we define functions h c 

' __4>{x) 
\ 

h aj (x) :- 





) 

- X 


0( aj ) - 


0(x) 


y/X - 
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I and hfy : (a 3 ,a j+1 ) 
fli-x < x < aj, 

aj < x < f3j, 



by the 



and 



(53) 



{ 0(x)~9(Pj) 
\[Pi - x 

\A - Pj ' 



ay < x < Pj, 



Pi < x < a 



3+1- 



Under the assumption of Condition [o] the definition ( 52 ) extends by continuity to x = ay and the definition 
( 53 1 extends by continuity to x = Pj ; moreover, these functions will all have one Lipschitz continuous 
derivative. Moreover, if x is bounded away from the support interval endpoints, h aj (x) and h aj (x) will have 
a second derivative that is also Lipschitz. This is shown in the Appendix in the case of G — but the 
same reasoning also works for G > 0. Note that the nonnegativity of the equilibrium measure implies that 
h aj {x) > and hpAx) > for ay < x < Pj, and the variational inequality (451 implies that h a Ax) < for 
Pj-i < x < aj and that hp, (x) < for Pj < x < ay+i. 

Now we may state the rest of the conditions that we impose on the external field V. 



Condition 2 (Strict inequalities and behavior at endpoints). For j 

for aj < x < Pj and that the functions h a 
(53 I satisfy the strict inequalities 



, G+ 1, we suppose that ip{x) > 
R defined by ( 52 ) and 



(54) 
and 
(55) 



h a . (x) < for Pj. 



< x < a 



j ■ 



hp j (x) > for aj < x < Pj 



(Pj-x,Pj) -> E and hp } : (aj,a j+1 ) - 

h a . (x) > for aj < x < Pj, h a . (aj) > 0. 

hp. (x) < for Pj <x < a j+1 , h' p (Pj) < 0. 



Condition 3 (Single interval of support w.l.o.g.). We assume that G = 0. 

The analysis for the case of G > (i.e. more than one interval comprising the support of fj,*) may be deduced 
in a straightforward manner from the case of G = (i.e. one interval comprising the support of /x*). So, 
in the course of our presentation of the details of the asymptotic analysis of the orthogonal polynomials, we 
will assume, without loss of generality, that G — and hence the equilibrium measure is supported on 
the single interval [a,P] = [ai,fii]. 

1.5. Statement of results. Because of the complex-conjugation symmetry p n (z*) — p n (z)* , we only need 
to present asymptotic formulae for the orthogonal polynomials in the upper half-plane. While our methods 
yield asymptotic formulae valid throughout the whole complex plane, we will restrict our attention to the 
regions £1+ and C+ n Sp as indicated in Figure [2] We focus on these regions for simplicity and also because 
these are most important for applications to random matrix theory. 




Figure 2. The regions f2± of the complex plane surround the interval (a,P). The square 
regions S a centered at a and Sp centered at P are shown with dashed boundaries. These 
squares have sides of length 26 for some 6 > 0, and the regions f2± each have vertical 
thickness 5. 

Theorem 1. Suppose the external field V satisfies Condition^through Condition^as described earlier. Let 
n, N — > oo so that N/n — > c with < c < oo. Then, with (a,P) representing the support of the equilibrium 
measure /i*, the following asymptotic descriptions are valid. 
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The asymptotic behavior of the leading coefficients K n ~i,n-i and K n ^ n is given by (2071. 

For z within f2 + but outside the squares S a and Sg, the orthogonal polynomials p n (z) and p n -\{z) 



possess the asymptotic description (212 1 and (213 1 respectively (see also (215 1 -(216) and (2211 

Pi;. 

For z within Sf)C\C+, the orthogonal polynomials p n (z) and p n ^\{z) possess the asymptotic descrip 



Hons (p25|-(p26| respectively (see also (p33|-p34| ; p37| -(p38| ; and ([241]) (p42|) y>. 

4. Asymptotic formulae for the derivatives of the orthogonal polynomials may be obtained by differen- 
tiating the leading-order asymptotics for the polynomials themselves, as described in Section^ 

These formulae describe the orthogonal polynomials in terms of the first column of the matrix A(z) given 
in (p59| , and the error terms therein are expressed in terms of the quantity A n defined by ( 202 1 . To mediate 



between the orthogonal polynomials contained in the first column of A(z) and the orthonormal polynomials 
Pn-i( z ) and p n (z), one must normalize by the leading coefficients K n —\ „_i and n nn respectively, whose 



asymptotic behavior for large n is given by (2071. As mentioned at the end of the previous subsection, the 



assumption that the support is a single interval is for convenience of presentation only. Theorem [T] may be 
easily extended to more general settings. As an example, it straightforward to carry out all the details if 
one assumes only Conditions [Oj [T] and [2] of the previous subsection. The following Theorems, describing the 
application of our results to random matrix theory, emphasize this point. 



Theorem 2. Suppose that the external field V satisfies Conditions'^ of subsection 1.4 Then Asymptotic 
Result ^from Section \l.l\ holds true. 



Proof. The calculations to deduce ( 20 1 from the asymptotic results concerning the orthogonal polynomials 
are by now standard, and we refer the reader to [16] for the details. (See also [6].) □ 



Theorem 3. Suppose that the external field V satisfies Conditions'^ of subsection 1.4 Then Asymptotic 
Result^from Section \l.l\ holds true. 

Proof. We again refer the reader to [TB] for the details of this calculation. (See also [5].) □ 

1.6. Notation. We will indicate complex conjugation with an asterisk: z* . All matrices are written boldface 
with the notable exception of the identity matrix I and the Pauli matrices: 



(56) 



0-1 := 



1 

1 



02 ■= 



-i 

1 



&3 := 



1 

-1 
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2. Orthogonal Polynomials, Riemann-Hilbert Problems, and Equilibrium Measures 

2.1. Characterization of orthogonal polynomials via a Riemann-Hilbert problem. Let N > be a 

parameter, and let V(x) be a real-valued function satisfying merely the conditions set down in the beginning 
of Section Q] 

The following Riemann-Hilbert problem |10j is known to characterize the polynomials {p n }^Lo orthonor- 
mal with respect to the measure vm given in ([!]), and defined by the conditions ([2|-(|3]). 

Riemann-Hilbert Problem 1. Find a 2 x 2 matrix A(z) = A(z; n,N) with the properties: 

Analyticity. A(z) is analytic for z € C \ R, and takes continuous boundary values A + (x), A_(a;) 
as z tends to x with i£l and z G C+, z G C_. 

Jump Condition. The boundary values are connected by the relation 
(57) A+(x) = A_(x) 
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I e -NV{x) 
1 



Normalization. The matrix A(z) is normalized at z = oo as follows: 
(58) fooA(*)(*~ n °„)=I. 

It was discovered in |10j that Riemann-Hilbert Problem [I] characterizes polynomials orthogonal with 
respect to vjq. The connection between these orthogonal polynomials and the solution of Riemann-Hilbert 
Problem [T] is the following: 



(59) 



A(z) 



V 



1 



-Pn(z) 



-2mK n - ltn -xP n -l(z) -K„_i,„_i 



Pn(s)e- NV ^ , 
as 

p n _ 1 (s)e- N ns) 



ds 



Note that ( 59 1 implies in particular that 

(6°) «n-l,n-l : 



lim z ( " 1) A 2 i(z) and «£, 

27TZ Z-+00 



— lim z- {n+1 ^ A^z)- 1 . 
2iri z— *oo 



These relationships provide a useful avenue for asymptotic analysis of the orthogonal polynomials in the 
limit n — > oo; it is sufficient to carry out a rigorous asymptotic analysis of Riemann-Hilbert Problem [T| 



2.2. Use of the equilibrium measure. Given the function g{z) defined by (46), we introduce an explicit 
change of dependent variable into Riemann-Hilbert Problem [T| Set 



(61) 



B(z) := e- rUa ^ 2 A(z)e- n9 ^ a3 e nia3 ^. 



It follows from the properties of A(z) set out in Riemann-Hilbert Problem [T] that the matrix B(z) is analytic 
for z £ C \ K and satisfies the normalization condition 



(62) 



lim B(z) 



The boundary values B+(x) and B_(x), taken on the real axis as z — ► x from the upper and lower half-planes 
respectively, are continuous functions of x S R related by the following jump condition: 



(63) 



= B_(x) 



-n(g + (x)-g-(x)) p n{g + {x)+g-{x)-cV{x)-t) 







c ,n{g+(x)-g-(x)) 







^in8(x) 



where 6{x) is defined by (49) and 4>{x) is defined by (51 1. From (48 1 and (49), we see that for x < a and 
x > [3 (recall that without loss of generality we are assuming that supp(?/>) = [a, 0\) this jump condition can 
be equivalently written in the form 



(64) 



B+(x) = B-(x) I ^ ), x < a or x > p. 



Similarly, from (44 1, we see that for a < x < (3 the jump condition takes the form 

!(x)j . a<x < (3. 



(65) 



, p — inO{x) 

B + (x)=B_(x) ( Q 



3. Extensions of 6{x) and <j>(x) 

In this section we will define extensions from certain intervals of the real axis of the functions 0(x) and 
(x). We shall assume throughout the conditions on the external field V described in the Introduction. 
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3.1. Existence of extensions with required properties. 

Lemma 1 (Extension of 0(x)). Suppose that c > is held fixed as n — > oo so that the real-valued function 
9(x) is independent of n. There exists a function Q(x,y) that satisfies the following: 

Property 1: Domain, smoothness, and boundary behavior. The function Q(x,y) is defined for 
a < x < (3 and \y\ < 8 for some 8 > 0. In its domain of definition, Q(x,y) and the partial derivatives 
Q x (x,y) and Q y (x,y) are all continuous and uniformly bounded. Moreover, if x + iy is bounded away from 
both a and (3, the second partial derivatives Q xx (x,y), Q xy (x,y), and Q yy (x,y) are also continuous and 
bounded. The function Q(x,y) is an extension of the real-valued function 0(x) in the sense that 

(66) Q(x,0) = 0{x), a<x<[3. 

Property 2: Behavior near the real axis. There exist finite constants K > and k > 0, such that the 
following three estimates hold true: 

(67) \dG(x,y)\ <K\y\\x + iy-a\ 1/2 \x + iy-[3\ 1/2 , a < x < (3, \y\ < 6. 

(68) Im(e(a;, y)) < -ky 3/2 , a < x < (3, < y < 6. 

(69) hn(Q(x,y)) > k\y\ 3/2 , a < x < (3, -5<y<0. 
Property 3: Behavior near a and (3. The function 

2n-Q(x,y) 



(70) G a (x,y) := , 

[x + iy — a) A i z 

extends continuously to x + iy = a with the limiting value G a (a, 0) = h' a (a) > 0. Moreover, 

(71) G a (x,y) = h' a (a) + 0(\x + iy-a\), asx + iy—>a, 
and 

(72) G a (x, ±(x - a)) = h' a (a), a<x<a + S. 
Similarly, the function 

(73) G,(x,y):= 



{p-{x + iy)fn 

extends continuously to x + iy = [3, and the limiting value Gp([3, 0) = — h'p(f3) > 0. Moreover, 

(74) G {x,y) = -h'^(3)+O{\x + iy-[3\), asx + iy^/3, 
and 

(75) Gf,(x,±(0 - x)) = -tip(0), (3-5<x<(3. 

Lemma 2 (Extension of 4>(x)). Suppose that c > is held fixed as n — > oo so that the real-valued function 
4>{x) is independent of n. Then there exists a function &(x,y) that satisfies the following: 
Property 1: Domain, smoothness, and boundary behavior. The function Q(x,y) is defined in two 
rectangles: R a given by a — 25 < x < a with < y < S and Rp given by (3 < x < (3 + 26 with —5<y<0 
for some S > 0. In its domain of definition, <&(x,y) and the partial derivatives <& x (x,y) and <& y (x,y) are all 
continuous and uniformly bounded. Moreover, ifx + iy is bounded away from both a and (3, the second partial 
derivatives $ xx (x,y), & xy (x,y), and & yy (x 7 y) are also continuous and bounded. The function &(x,y) is an 
extension of the real-valued function 4>{x, y) in the sense that 

(76) $(x,0) = <j>(x), (3 < x < (3 + 28 and a -28 < x < a. 

Property 2: Behavior near the real axis. There exist finite constants K > and k > 0, such that the 
following estimates hold true: 

(77) \d${x,y)\ <K\y\\x + iy-a\ 1/2 \x + iy- f3\ 1/2 , (x, y) G R a U R , 

(78) Re(®(x,y)) > k\x + iy - a\ 3/2 for (x,y) e R a and Re($(x, y)) > k\x + iy - /?| 3/2 for (x, y) € Rp. 
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(Note that from (78 1, the weaker inequality 

(79) Re(<I>(x,y)) >fc|y| 3 / 2 , 

follows immediately.) 

Property 3: Behavior near a and (3. The function 



(x,y) G R a U Rp 



(80) 



H a (x,y) 



(a — (x + iy)) 3 / 2 

extends continuously to x + iy = a with the limiting value H a (a, 0) = h' a (ct) > 0. Moreover, 



(81) 
and 
(82) 

Similarly, the function 
(83) 



#aOEj y) = h' a (a) + 0(\x + iy - a\), as x + iy 
H a (x, ±(a — x)) = h' a (a), a — 5<x<a. 

Hp{x,y) — 



((x + iy)- 0) 3 / 2 

extends continuously to x + iy = (3, and the limiting value satisfies Hp(f3, 0) = —h'a((3) > 0. Moreover, 

(84) H (x, y) = -h'pifi) + 0(\x + iy - /3| fe - 2 ), as x + iy - (3, 
and 

(85) Hp(x,±(x-(3)) = -tip(/3), (3<x<(3 + 8. 

3.2. Proofs of Lemmas [l] and [2} In this subsection we construct suitable extensions Q(x,y) and &(x,y) 
by further developing a strategy used in |15j . We will use the following notation generalizing the "bump" 
function B(-) introduced in Section [I] for an interval [a, b], 



(86) 



B[ a ,b](x) := B 



x — a 
b — a 



This C(°°)(R) function maps M onto the interval [0,1] with B^ ab ^(x) = for x < a and Br a u(x) = 1 for 
x > b. 

3.2.1. Proof of Lemma^ extension of 9 (x). 

Definition of the extension 0(x,y). First, we define functions & a ,o( x ,y) an d ®p,o(x,y) by 

(87) <da y0 (x, y) := 2tt - ((x + iy) - a) 1 !" 1 [h a (x) + i(h a (x + y) - h a (x))] , x < [3 and x + y < (3, 
and 

6/3,0(2:, y) := ((3 — (x + iy)) 1 / 2 [hp(x) + i(hp(x + y) — hp(x))] , x > a and x + y > a. 



Here, the function h a (x) = h ai (x) is defined by (52 1 for —00 < x < (3 and the function hp(x) — hp 1 (x) is 
defined by (53 1 for a < x < +00. 

< Remark: The extensions of the functions h a and hp within the square brackets in (87 1 and (88 1 
respectively are Cartesian versions of the polar-coordinate extensions discussed in |15j , further generalized 
with the use of difference quotients in place of derivatives. > 

Next, we define analytic approximations of Q ai o(x, y) and Qp t o(x, y) valid for x + iy « a and x + iy ~ (3 
respectively: 

(89) ®l°o(x,y) :=2ir-ti a (a)((x + iy)-a) 3/2 , x > a or y ? 0, 
and 

(90) e^(x,y):=-hp(f3)(P-(x + iy)) 3 / 2 , x<(3ory^0. 

In precisely the spirit of the simple example described in Section [T] we may combine the two types of 
extensions with the help of an appropriate angular bump function: 



(91) 



Q a (x,y) := B 



14 



1 - B 



® a ,o(x,y), 



and 



(92) Q (x,y):=B 
For short we will occasionally write 



y 

x- /3 



e£o(*,w) + 



1 - B 



y 

x- 



0/3,o(x,y). 



(93) 

Finally, letting 
(94) 



,111:4 



a '■— B 



and B 



ang,/3 



:= B 



2-/3 



and 



6:=/3--(/3-a) 



so that [a, b] C (a,/3), we may smoothly glue these two extensions together through the vertical strip 
a < x < b in the {x, y)-plane: 



(95) 



®(x,y) := B [a ^ b] (x)Qp(x,y) + [l - B [aM (x)] <d a {x,y). 



This will be our extension of the function 0(x) from the interior (a, (3) of the support interval. Taking into 
account the supports of Br a w (x) and 1 — B^ a ^ (x) and comparing with the regions of definition of Q at o(x, y), 
Qf3fi(x,y), Qa°l(x,y), and 6j3°o(x, y), we see that whenever S < ((3 — a)/3, &(x,y) is well-defined on the 
rectangle R given by the inequalities a < x < (3 and \y\ < 5. 

< Remark: It turns out that if we replace both angular bump functions -Bang, a and i? a ng,/3 by the 
constant function 5=1, then the extension obtained only involves the functions Q a fi{x,y) and Q[j$(x,y) 
glued together through the vertical strip a < x < b, and this simpler function satisfies all of the desired 



properties with the exception of (75 1 and (72 1 from Property 3. The purpose of the angular bump functions 



is to smoothly deform the simpler extension into one that satisfies these additional conditions (without 
ruining any of the other conditions, of course). > 

Confirmation of Property 1. To confirm Property 1, we note that under the assumptions in force, both func- 
tions h a and hp have one Lipschitz continuous derivative throughout their respective domains of definition, 
which implies that d x <d at o(x, y), d y Q ai o(x,y), d x Qpfi{x,y), and dyQp^ (x, y) are all continuous and uni- 
formly bounded throughout the rectangle R of definition of Q(x,y). Since 6^°q(x, y) and 6^°q(x, y) are 
analytic functions for (x, y) € R, their first partial derivatives are certainly continuous. Then, since Q(x,y) 
is constructed from these more elementary functions with the help of C' 00 -' bump functions, it is clear that 
Q(x,y), Q x (x,y), and Q y (x,y) are all continuous an uniformly bounded throughout R. As h a {x) and hp(x) 
have a second Lipschitz derivative for x bounded away from a and (3, similar arguments show that Q xx (x, y), 
®xy(%i y)i and Q yy {x, y) are continuous and bounded for x + iy bounded away from a and (3. Furthermore, 
for a < x < (3, 



(96) 



Q(x, 0) = B [aM (x)e l3 (x 1 0) + [1 - B [aM {x)] G a (x, 0) 

= B [aM (x)e 0fi {x, 0) + [1 - B [aM (x)] Q a , {x, 0) 
= B [aM {x)6{x) + [l-B [aM (x)] 9(x) 
= 0(x), 



so Q(x,y) is indeed an extension of 0(x) from the interval (a, 13) to the rectangle R. 
Confirmation of Property 2. To confirm Property 2, first note that from (951 we have 



(97) 



dQ = dB [a>b] ■ (Op - Q a ) + B [a . b] dQp + [1 - B [aM ] de a . 



Now, since & a (x,y) and Qp{x,y) are both extensions from (a,/3) of the same function 9(x), and since they 
are both uniformly Lipschitz for x + iy bounded away from a and (3 (this is where dBt a ^ is nonzero), the 
first term on the right-hand side is supported in a < x < b and is C(|y|). Therefore we certainly have 



(98) 



\d B [a,b] ■ (6/3 



e a )\<K\y\\x- 
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iy- 



al^lx- 



1/2 



for some constant K > 0. It therefore remains to estimate dQ a for x + iy bounded away from (3 and 30/3 
for x + iy bounded away from a. Since dQ^^x, y) = and dO^fa, y) = 0, we see that 

dQ a = dB anga ■ (6^°q — 6 Qj0 ) + [1 — -S ang , Q ] dQ afi 



(99) 



Now, by direct calculation 



96,3 = DB^.p ■ (6*° - 6/3,0) + [1 - B WSil3 ] de 



0,0- 



(100) 



dB ; 



ang,/3 



B' 



y 



d 



V 



B' 



!J 



1 \x + iy — (3\ 

2 (x-/3) 2 ' 



x — (3 J x — P \ x — P 

and since the inequality |a; — /3| > \y\ holds whereever the derivative of the bump function in this formula is 
nonzero, 



(101) 



\dB. 



sag, /3 



< 



B' 



y 



x-P 



\x + iy- P\ 
(x - P) 2 + y 2 





V 


) 




x-P 





1 



\x + iy- P\ 



Also, since for any w, 

(102) h f3 (w)= I ti ( s )d S = h'p(P)(w-P)+ I [hfa) - tipffl] ds, 











we have 

(103) e%{x,y)-e fii0 {x,y) 

= -(P-(x + zy)) 1 / 2 



(1 - i) / - h'pifl] ds + i [h' {s) - h'piP)] ds 



x+y 



so, since hp(x) has one Lipschitz continuous derivative, there are constants K\ > and K 2 > such that 

\e h p °\x, y) - 6/3,0(2;, y)\<\x-P + iy\ x ' 2 [K^x - P) 2 + K 2 (x -p + y) 2 } 

(104) =\ x -p + iy \^ 2 [2K(x - Pf + 2K 2 y(x - p) + K 2 y 2 ] 

< \x-P + iy\^ 2 [2K(x~P) 2 + 2K 2 \y\\x-P\ + K 2 y 2 ] , 

where K := (K\ + K 2 ) /2. Using again the inequality \x — P\ > \y\ (since we are going to multiply by dB ang ^), 
we therefore have 



(105) 



\G h pf(x, y) - 6/3,0(2;, y)\ < \x-P + iy\ x ' 2 [2K(x - P) 2 + 3K 2 \y\ \x - p\] . 



Since we have both \B'(t)\ < C and \B'(t)\ < C\t\ for some constant C > 0, it follows from ( flOlj ) and ( |T05 l 
that 



(106) 



as a „g,/3 • (e^°o - 6/3,0) = O (\x + iy - pr 1/2 \y\\x - p\) 

= o(\y\\x + iy-p\ 1 l 2 ) , 

where we have used \x — P\ < \x + iy — P\ in the last step. Furthermore, 
(107) 

dGp <0 (x, y) = (p-ix+ty)) 1 ^ [hp(x) + i(hp(x + y) - hp(x))] = l -{i-l){P-{x+ l y)) 1 ' 2 %{x + y) - h' (x)] , 
so since h'p(x) is uniformly Lipschitz near P and 1 — -B ang ,/3 is bounded, we also have 

(108) [1 - B ang ,/3] 96,3,0 - O (\y\\x + iy- p\ 1 ' 2 ) . 

Therefore, for x + iy bounded away from a we have 



(109) 



\dGp(x, y) I < K\y\ \x + iy- a\ 1/2 \x + iy- P\ 1/2 



for some constant K > 0. In a completely analogous fashion we see that for x + iy bounded away from P we 
have 



(110) 



\dQ a (x,y)\ < K\y\\x + iy- a\ l ' 2 \x + iy - p] 1 ' 2 . 
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Combining these results with ( 98 1 we complete the proof that 

(111) \dG(x,y)\ < K\y\\x + iy - a\ 1/2 \x + iy - f3\ 1/2 . 

Now consider lm(0(a;, y)). Since all of the bump functions Blm, -Bang, a, and B aTL g,p are real-valued, it 
will suffice to analyze Im^^o^) y)) and Im(6^°Q(a;, y)) for x + iy e R bounded away from (3 and to analyze 
lva{Qpfl{x,y)) and Im^^^x, y)) for x+iy £ R bounded away from a. Writing /3— (x+iy) = \(3— (x+iy)\e 1 ^ 
with |0| < 7r/2, we have the exact formulae 



(112) Im(e Pfi (x,y)) = \l3-(x + iy)\ 1 / 2 y 

hp{x + y) - hp(x) 



<A _ 1 hp(x) 
2 J 2 ' P-x 



1 - tan^ 



and 
(113) 



Ira(Q h p %x,y)) = \p-(x + l y)\ 1 ^y 



sec l- 



Since Condition [2] requires that h!p{f3) < 0, the condition that \cf>\ < n/2 immediately implies that 



(114) 



i^(/J) (4008 (|) 



seel- 



< \h' p {P) < 0. 



To analyze the terms in the square brackets in (112) requires a little more work. Suppose first that \/3— (x- 
iy)\ < e\. Then, as ei — > 0, we have both 



(115) 



M* + y)-M«) _ w<0 and _i.^M^i w<0> 



where the inequalities follow from Condition[2] So ei may be taken to be small enough that \(3— (x + iy)\ < e\ 
implies both 



(116) 



h fJ (x + y)-h fj (x) < 1, ^ 
y 2 



1 . < 1 W < 0. 

2 p — x A 1 



Therefore, since \<j>\ < tt/2 implies both 

(H7) ^ 

we see that \0 — {x + iy)\ < t\ implies 
hp(x + y) - hp{x) 



2 <C ° S V2 



and < 1 — tan 



(118) 



V 



h (x) f 
2J 2 j3- x \ 



tan" 



V2 U 

4 



On the other hand, if we suppose that \/3 — (x + iy)\ > e±/2, but that x > a and \y\ < €2 for some £2 > 0, 
then as £2 — > we have both 



(119) 



and W X + y) -W x) -+hfr). 



Now, from ( 53 1 , for a<a<ir</3we have 

1 hp(x) 



(120) 



h'Jx) 



y 



V{x) 



2irip(x) 



2(3-x y/W^x yfW^x 



< -h 



with some constant k\ > as a consequence of Condition [2] and the square-root vanishing of ip(x) as x | P- 
Therefore, by choosing £2 sufficiently small we will have 



as long as \/3 — (x + iy)\ > e\/2, x > a, and \y\ < To combine these estimates, note that if 6 > is 
sufficiently small, the part of the rectangle R given by the inequalities a < a < x < (3 and \y\ < 6 consists 
of points (x, y) for which either \(3 — (x + iy) \ < t\ or |y| < e 2 , so ( 118 1 and ( 121 1 may be combined to give 



(122) 



hp(x + y) - hp(x) 



for (x,y) € R, where 
(123) 




tan 



4 ,^03)1,2*1 



< -k 2 < 



Finally, we may combine ( 114[ ) with ( 122 1 to find that for (x, y) € R with x > a, 



iy-P\ 



1/2| 
|V2| 



Im(6 /3i o(a;,2/)) < -fay\x 
Im(e^(x,y))<-k 3 y\x- 

lm(Q,3 y0 (x,y)) > -k 3 y\x- 
Im(6^(x, y)) > -k 3 y\x + iy - 0\V 2 \x 



-iy-PV 
^y-P\ 1/2 \ 



- w 

- iy 

iy - 
iy - 



a 



(124) 
and 
(125) 
where 

(126) , ,. _ , _ . 

y/(/3 - a) 2 + S 2 

Completely analogous arguments show that for (x,y) S R with x < b, 



1 1/2 
1 1/2 



1/2 
1 1/2 



if V > 



if y < 0, 



(127) 

and 

(128) 



Im(6 ai0 (x, y)) < -k 4 y\x + iy - /3\ 1/2 \x + iy - 
Im(e^(z, y)) < -k 4 y\x + iy - p^x + iy - 

Im(6 Qi0 (a;, y)) > ~k 4 y\x + iy - /3\ 1/2 \x + iy - 
Im(e h °l(x,y)) > ~k 4 y\x + iy - (3\ x ' 2 \x + iy - 



1 1/2 
1 1/2 



1 1/2 
1 1/2 



if 2/ > 



if y < 0, 



for some constant fc 4 > 0. Letting fc = mhi{k 3 , k 4 } > and using the fact that <d(x,y) is a convex 
combination of Im(6 Q! o(a;, y)), Im(6^°Q(a;, y)), Im(®p t o(x,y)), and Im^^^x, y)) through the various bump 
functions involved in the definition, it follows that 



(129) 



Im(6(a;, y)) < -ky\x + iy - (3\ 1/2 \x + iy - a\ 1/2 , y>0 
Im(6(x, y)) > -ky\x + iy~ f3\ 1/2 \x + iy~ a\ 1/2 , y < 



holds for (x, y) € R if the thickness parameter S of the rectangle R is sufficiently small. Now, if x < (a + /3)/2 
then \x + iy ~~ (3\ is bounded away from zero while \x + iy — a\ > \y\, and if x > (a + (3) /2 then \x + iy — a\ is 
bounded away from zero while \x + iy — (3\ > \y\. Combining these observations with ( |129 1 yields (68 1 and 
([69]). 

Confirmation of Property 3. To confirm that @(x,y) satisfies Property 3, note that 

@/3,o0e, y) _ hp(x) + i{hp{x + y)- h (x)) 



(130) 



Gp,o{x,y) 



((3 - (x + iy)) 3 / 2 
is continuous near (x, y) = (/?, 0) and satisfies 
(131) Gp, (x,y) = -ti l3 (P)+O(\l3- 

because h^x) is Lipschitz continuous. Similarly, 



(132) 



09 - (x + iy)f/ 2 
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(3 - (x + iy) 
(x + iy)\) 

= -h'pifl) > 



so since near x + iy = f3 (that is, for x > b) Gp(x,y) is a convex combination of Gpfl(x,y) and Go°q(:z;, y), 



the requirement (74 1 on Gp{x,y) given in Property 3 is met. And since for \y\ > (3 — x we have Gp{x,y) — 
G^°l(x,y), we also confirm the requirement (751. Similar calculations show that G a (x,y) satisfies the re- 
quirements fill and (172). 



3.2.2. Proof of Lemma []| - extension of 4>(x). The construction of a suitable extension of </>(x) follows the 
same general procedure as the construction above of <d(x,y). We give all details of the construction, after 
which it is straightforward to follow the reasoning given in the proof of Lemma [l] to establish that Properties 
1, 2, and 3 are satisfied. 
We first define 

(133) $ a fi(x, y) :— —(a — (x + iy)) 1 ^ 2 [h a (x) + i (h a (x + y) — h a {x))] , x < (3 and x + y < (3, 
and 

(134) $/9,o (a, y) ■= -((x + iy) - (3) 1/2 [hp{x) + i (h p (x + y) - hp(x))] , x > a and x + y > a. 
The analytic approximations of these functions valid for x + iy w a and x + iy « [3 respectively are 

(135) Os-f) ^h' a (a)(a~(x + iy)f/ 2 , x<aory^0, 
and 

(136) <S> h ° l o (x,y):=-h' (f3)((x + w)-f3) 3/2 , x>(3ory^0. 

Since the rectangles R a and Rp are disjoint, there is no need to merge functions defined near x + iy = a 
with functions defined near x + iy = (3, so we may simply define 



(137) $(x,y) 



D 



D 



IS 



x-(3 



1 - B 



1 - B 



V 



x- (3 



®afi{x,y), (x,y) e R a 
$0,0(2;, y), (x,y)€Rp. 



4. An equivalent Riemann-Hilbert-3 problem. 

The jump condition satisfied by the boundary values taken by B(z) on (a, /3) can be written in the factored 
form: 



(138) 



B + 0r)=B_(x) 



1 



Consider the contour E illustrated in Figure [3] 







■iB(x) ^ 



1 

-1 



1 

e -inO(x) ^ 





Figure 3. The oriented contour £ consists of the real intervals (—00, a — 26), (a,/3), and 
(/3 + 2<5, +00), along with the indicated contour segments S,, connecting a — 26 to a, £+ and 
£_ connecting a to f3, and E/3 connecting (3 to (3 + 26. All contour segments are oriented 
left-to-right, and all nonhorizontal segments have slopes ±1. 

Let Q(x,y) be any extension of 0(x) having all three properties described in Lemma [TJ and let <f>(x,y) 
be any extension of (j)( x ) having all three properties described in Lemma [2j We define a matrix D(x,y) for 
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z = x + iy G C\E relative to the domains f2_, f2 Q , and O^g shown as shaded regions in Figure [3] as 
follows. Set 

(139) Y>(x,y) :=B(x + iy)(^_^ i l e{x y) , x + iyett + , 

(140) T>(x,y):=B(x + iy)(^ ein Q {Xty) , x + iy€fl_, 

(141) D(a:,y) :=B(ir + iy) I Q x J, x-MyGft«, 

(142) D(x,j/) := B(x + iy) \ t x , x + iyeflp, 

and for all remaining z £ C \ E, we set D(x, y) := B(x + iy). 

Because it is explicitly related to B(x + iy) and hence to A.(x + iy), the matrix D(x, y) will solve Riemann- 
Hilbert-3 problem [l] to be defined below. Define the jump matrix Vd(z) for z€Eas follows: 

(143) Vd(«):=(q x ), z = x <a-25 and z = X > /3 + 25, 

e j ) » 2 = z + %e£ a u£/?, 

(145) V D (z):= ^ JV z = xG(a,/3), 

(146) V D (z) := r_i„i( X)1 ,) z = i + iy€S + , 

(147) V D (z) := f^e^.i,) ' Z = x + * y G 
Also, define the auxiliary matrix Wo(x, y) as follows. 

(148) Wo ^'^ := (me-» e ^)ae(x, y ) o) ' * + ^ eQ +> 

(149) W (x, y ):= (^^(^j J)' s + H^"-' 

(150) W (x,y) := I Q Q x + zt/Gft Q , 

(151) W (x,v):=l Q v ' y; ), x + zye^. 

For all remaining (x,y) G M 2 , we set W (x,y) :— 0. Note that W (x,y) so-defined is compactly supported. 
From the properties of the matrix B(z) inherited via the substitution (61 1 from properties of the matrix 
A(z) contained in the statement of Ricmann-Hilbert Problem]!] it follows that D(x,y) solves the following 
hybrid Ricmann-Hilbert-9 problem: 

Riemann-Hilbert-9 Problem 1. Find a 2 x 2 matrix T>(x,y) with the properties: 

Continuity. D(x, y) is a continuous function of x and y for x + iy G C \ £ taking continuous 
boundary values D + (x,y) (respectively D_(x, y)) on £ from the left (respectively right). 
Jump Conditions. The boundary values are connected by the relation 

(152) D+(x,y) = D_(x,y)V D (x + iy), x + iyeS. 
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Deviation From Analyticity. For x + iy € C, 

(153) dT>(x,y) =T>(x,y)W (x,y). 

(Note that in particular d~D(x, y) = for x + iy ^ f2 + U f2_ U fl a U O^.,) 
Normalization. The matrix T)(x,y) is normalized as follows: 

(154) lim T>(x,y) = L 

x,y — *oo 

5. Construction of a Global Approximation to T>(x,y) 

In this section we will build a global approximation to D(x, y) by considering a Ricmann-Hilbcrt problem 
obtained from Riemann-Hilbert-9 problem [I] by ignoring the "<9 component" of the problem: 

Riemann-Hilbert Problem 2. Find a 2 x 2 matrix D(z) urai/i i/ie properties: 

Analyticity. D(z) is an analytic function for zeC\S taking continuous boundary values D+(z) 

(respectively T)^(z)) on S /rom i/ie Ze/( (respectively right). 

Jump Conditions. The boundary values are connected by the relation 

(155) D+(z) = D_(z)V D (x,y), z = a: + iye£. 
Normalization. The matrix D(z) is normalized as follows: 

(156) lim D(z) = I. 

2 — *00 

Riemann-Hilbert Problem [2] has been obtained from Riemann-Hilbert-9 Problem [T] in an ad-hoc fashion, 
and even though T)(x,y) clearly exists, it is not immediately clear that a solution D(z) to Riemann-Hilbert 
Problcm|2]exists. Theorem [4] below asserts that a unique solution exists, and describes important asymptotic 
properties of D(z). The Theorem describes the asymptotic behavior of the solution in three different regions 
of the complex plane: two square domains S a and Sp of side-length 26, with S a centered at a and Sp 
centered at j3, and one exterior domain, C \ (S a U Sp). We further subdivide each square into four regions 
according to the contour E as indicated in Figure |4j 




Figure 4. The squares S a (left) and Sp (right), each subdivided into four regions as indicated. 

Note that according to Property 3 in Lemma [l] and Lemma [2j the restriction of the jump matrix ~Vr>( z ) 
to S R S a and S n Sp is piecewise analytic. Indeed, if we define 

(157) u a = u a (z) := [h' a (a)} 2/3 (a - z), and up = up{z) := [-h' p (P)] 2/3 (z - 0), 
with the positive two-thirds power meant in each case, then we have 

(158) V D (z) =r e 1 , zeS Q nS n , 
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H')'M V D c): JV - e m. nn.s;, 

(160) V D ( Z ) - ^^3/, JV zeE±n 5„, 

and 

(i -JlUa \ 

J j J. ^^ns 3 , 

(162) V D (z) = , z S (a,/?) n % 

(163) v D (z) = ^ , z e E± n fy. 

For all z e £ outside of the squares S a and S^, with the notable exception of the interval (a + 6, — 5) 
where Vd(z) is a constant matrix, the jump matrix Vrj(z) decays exponentially to the identity matrix as 
n — ► oo, as a consequence of both the variational inequality 4>{x) > for x < a — 26 and x > /3+ 26 and also 
the inequalities on lm(0(x,?/)) in Property 2 of Lemma [7J and the inequality on Re(<i>(a;, y)) in Property 2 
of Lemma [21 

In a way that is by now quite standard (see [6j[7]), these facts suggest an explicit model for D(z) that we 
will call D(z) and that we will now define. Let 7(2) be the function analytic in C \ [<x,0\ determined by the 
conditions 

(164) 7 (z) 4 = Z —^- , and lim 7 (z) = 1, 

z — a z— »oo 

and let U denote the unitary eigenvector matrix for Vd(z) on (a,j3): 

1 /e - "/ 4 e i,r / 4 

(165) U:= ^U W4 e—/ 4 
Then, we set 

K7W + 70O- 1 ) s(7W"7( Z )-T 

(166) D(z) :=U 7 (z) CT3 Ut = ( |, 2 gC\(S a US^, 

i( 7 (z)- 7 (z)- 1 ) KtW+tW" 1 ), 



2; 
CT3/6 



(167) 



D(z) := -V2ttU 3 7 (z)^ UQ (z)^/ 4 a 2 MK(z)) f 7 3 e"^( z ) 3/2CT3 / 2 , z e S a , 

and 

(168) D(z) := V2^U (~ \ * 7 (z) CT3 U/ 3(z)- ff3 / 4 M(u /3 (z))e™ u ' 3(z)3/2ff3/2 , z € fy, 

where M(u) is defined as follows with £ := (^p) 2 ^ 3 it. 

... , fe-^Ai'tt) e n "/i2 A i' (ee- 2 "/ 3 )\ tt 3tt 

(170) M(u) := y e - 7 ^/i 2 Ai /^ e 2^/3< e ^/ 12 Ai (Ce" 2 ™/ 3 ) J ' T < ai ' g(u) < 

/ e iiWi2 Ai ' (^-aws) e 7Wi2 A i'fe 2 "/3)\ 3^ 

(171) M(u) := e47r/12A . ^_ 2ri/3 ) e 5 -/ 12 Ai (£e 2W3 ) J ' ^ K ai ' g(u) < 

... , /e- 3 "/ 4 Ai' (0 e 7 "/ 12 Ai'(ee 2 ^ 3 )\ 3^ . . tt 

(172) M(«) := e _ 47r/4A . W e57ri/12A . j , __ < arg(.) < --. 



Here Ai(£) denotes the Airy function, the unique solution of y" — £y = with the asymptotic behavior 

(173) Ai(0 = ^ i7i ^(l + O(|^|- 3 / 2 )) and Ai'(fl = (1 + 0(1^)) 

as £ — > oo with arg(£) € (— 7r, 7r). 

< Remark: For those readers familiar with the notation of the paper [7] we make the following clarifica- 
tion. The matrix M(u) defined here may be expressed in the form 

(174) M ^=i(? V)p((!) 2/3 ^^" 3/2ct3/2 



where the local parametrix P(£) is as defined in [7], equations (1.36)-(1.40). > 
The point of introducing the matrix D(z) is the following. 

Theorem 4. Assume the conditions on the external field V stated in the Introduction. Let n, N — ► oo so that 
N/n — * c lozi/i < c < oo. TTien for n sufficiently large, there is a unique solution D(z) to Riemann-Hilbert 
Problem^ which possesses the following global asymptotic description: 




(175) D(z) = 

uniformly with respect (ozGC as n — > oo . 

Proof. Let F(z) := D(2)D(z)^ 1 . It is easy to see from the properties of D(z) required by the conditions 
of Riemann-Hilbert Problem^ and the explicit formulae given for the matrix D(z) in various parts of the 
complex plane that F(z) is a matrix that is required to have the following properties. Firstly, F(z) is analytic 
at least for z £ C \ E F , where E F is the union of E and the boundaries of the square regions S a and Sp, and 
F(z) takes continuous boundary values on E F . Secondly, the boundary values satisfy F + (z) = F_(z)Vf(z) 
for some jump matrix function Vf(z) defined on E F that is explicitly calculable in terms of Vd(z) and the 
boundary values taken on E F by F)(z). Thirdly, F(z) must tend to the identity matrix as z — > oo. In other 
words, these three facts show that F(z) satisfies its own Riemann-Hilbert problem. 

The Riemann-Hilbert problem satisfied by F(z) is of a particularly convenient type: it is a "small-norm" 
problem in the sense that the jump matrix Vp(z) is a small perturbation of the identity matrix in a suitable 
space of matrix-valued functions on the contour E F . In fact, it is easy to check by direct calculation that 
Vf(z) = I for E F n (S a U Sp). This is a direct consequence of Property 3 in Lemma [l] and Lemma [2] 
characterizing respectively the extensions 0(x, y) and <j>(x,y) on these portions of the contour E, and of the 
identity 

(176) Ai(f ) + e - 27ri/3 Ai(£e- 27 "/ 3 ) + e 2 "/ 3 Ai(£e 27ri/3 ) = 0. 

An even easier calculation shows that ~Vp(z) = I for a + 5<z<(3 — 6. With the use of the asymptotic 
formulae (]173[), one sees that on the boundaries of the two squares S a and Sp, Vp(z)-I is uniformly 0(n _1 ), 



and on all remaining parts of E F one finds (in part by the estimates ( 68 1 — ( 69 1 on lm(Q(x,y)) in Property 2 
of Lemma [TJ and the estimate (78 1 on Re(Q(x,y)) in Property 2 of Lemma |2| that Vp(z) — I is uniformly 



exponentially small as n — > 00 and also decays rapidly as z — > 00. 

Since for our purposes we need to control the size of F(z) — I right up to the contour E F , we need to 
formulate the Riemann-Hilbert problem for F(z) in an appropriate space in which the boundary values 
of F(z) are Holder continuous with some exponent a £ (0,1]. To do this we need to observe that as a 
consequence of our assumptions on the external field V and the corresponding smoothness of Q(x,y) and 
<I>(x, y) described in Property 1 of Lemma [TJ and Lemma [2j and also as a consequence of the piecewise 
analyticity of the comparison matrix D(z), the jump matrix Vp(z) is sufficiently smooth on a sufficiently 
(piecewise) smooth contour that the Holder version of the small-norm theory applies. The result is that as 
n — > 00, F(z) exists uniquely in the space of matrices with Holder-continuous boundary values, and also 
F(z) — I = C(n _1 ) holds uniformly throughout the complex z-plane. 

A detailed account of the existence theory for small-norm Riemann-Hilbert problems is discussed, for 
example, in [S]. Specific information relevant to the application of small-norm theory in Holder spaces can 
be found in Appendix A of [11] . □ 
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An important property of T)(z) is that for all z where it is defined, det(D(z)) = 1. Therefore, D(z) and 
its inverse D(z) _1 have comparable bounds in any matrix norm. From ( |166 1 one may then see that D(z) 
and its inverse are bounded as n — > oo uniformly for z G C \ (S a U Sp). On the other hand, from (168) and 
([167} together with the definition (fl69|-([T72| of M(u) one sees that D(z) = £>(n 1/6 ) and D(z)" 1 = 0{n 1 / & ) 



hold for z G S a U Sp, although for our purposes a more useful estimate coming from the same formulae is 
that 

(177) D(z) = 0(|z-a|- 1 / 4 |z-/3|- 1 / 4 ) and D(z)- 1 = <D(\z - a\- 1/4 \z - /3\~ 1/4 ) 

holds uniformly for z in bounded sets (the constants implicit in the order relations are independent of both 
n and z). 

The construction of D(z) is one part of the argument where the details are somewhat different for G > 
than for G = 0. To handle the case with more than one interval of support one must replace the definition of 
D(z) for z G C \ (S a U Sp) with a matrix constructed from Riemann theta functions for hyperelliptic curves 
of nonzero genus modeled by two copies of the complex z-plane identified along cuts made on the real axis in 
the support intervals of the equilibrium measure /i*. Full details may be found, for example, in [TJ. The key 
property of uniform boundedness of D(z) away from the endpoints of the support intervals remains valid in 
this more general case. 

6. A d Problem and Existence Theorem 
Having constructed D(x + iy), we now define E(x,y) via 

(178) E(x,y):=T)(x,y)t>(x + iy)- 1 . 

It is immediately clear that E(x, y) is continuous in C. It is straightforward to compute the d derivative of 
E(x,y), and we learn that E(x,y) solves d Problem [T] below. We define a "dressed" version of the matrix 
Wj(i, y) as follows: 

(179) W(s,y) :=t)(x + iy)W (x,y)t)(x + iy)- 1 , (x,y)eR 2 . 

d Problem 1. Find a 2 x 2 matrix E(x, y) with the properties: 

Continuity. E(x, y) is a continuous function of x and y for (x,y) G K 2 . 
Deviation From Analyticity. For (x,y) G M 2 , 

(180) dE(x,y)=E(x,y)W(x,y). 

(Note that in particular <9E(x, y) = for (x, y) ^ f2+ U U fl a U Sip.) 
Normalization. The matrix E(x,y) is normalized as follows: 



(181) 



lim E(x,y)=I. 

x,y — >oo 



In view of the normalization condition ( 181 ) and the fact that W(x, y) = outside some compact set, we 

Cauchy k 

/CE(x,y), 



may invert the d operator in ( 180 1 with the help of the Cauchy kernel 
(182) 
where 

1 

7T 



(183) 



JCE(x,y) :-- 



E(x,y) 



((u + iv) — (x + iy)) 1 E(m, u)W(it, v) du dv. 



It is a basic fact that if E(x,y) satisfies the integral equation ( |182[ ) then E(x,y) also solves d Problem [lj in 
fact the integral equation ( 182| is equivalent to d Problem [T| 

In this section we will show that the integral operator /C, when considered in the space L°°(IR 2 ), has norm 
bounded by Cn _1 / 3 log(n) for some C > 0. This implies that the integral equation (1821 may be solved by 
Neumann series. 



The strategy to prove this is quite straightforward: because the singularity of the Cauchy kernel is 
integrable in K 2 , the basic estimate is: 

1 J[ ll w K u ) 
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(184) 



^H|| L oo( R 2) < ||H|| L =c( K 2) 



sup 

(x, y )em 



\((u + iv) - (x + iy))\ 



du dv 



where ||W(u, u)|| is a pointwise matrix norm, i.e. a norm of the matrix W evaluated at (u,v). Since W 
is uniformly bounded and has compact support, this immediately implies that K. is a bounded operator on 
L°°(M. 2 ). The goal is then to prove that the integral appearing on the right hand side of (184) is small, and 
for this the (u,v) dependence of ||W(u,w)|| will be essential. 

Theorem 5. There is a unique solution E to d Problem^ which possesses the following uniform asymptotic 
description, valid for all (x,y) £ M 2 : 

(185) e(*,») = i+V log(n) ^ 



1 / 3 y/l + x 2 + y 2 ) ' 

where the constant implicit in the order notation is independent of n and z and depends only on the external 
field V and the constant c. 



Proof. We begin by describing the asymptotic behavior of W(x, y). According to ( 179 1, W(x, y) is obtained 
from Wo(x, y) by conjugation, so we start by making the following two observations about Wo(x, y). Firstly, 

supp(Wo) = n + u u n a u n 

(186) ' , 

C B := {(x,y) e M 2 : a - 26 < x < /3 + 25, \y\ < 5}. 

Secondly, as a consequence of the definition of Wo(x,y) given in ( 148| — (151 1 and Property 2 in Lemma [l] 
and Lemma [2] describing d<d(x,y), Im(6(x,y)), d&(x,y), and Re($(x, y)) for (x,y) in the support of Wo, 
we may assume that 

(187) ||W (tt,t>)|| < Kne~ knM3/2 \v\\w-a\ 1/2 \w- ft 1 ? 2 , w = u + iveB, 

for some constants K > and k > 0. Using the definition ( 179 1 of W(x, y) in terms of Wo(x, y), Theorem|4] 
together with (1771 implies that 

(188) ||W(u,u)|| < Kn\v\e- kn ^' 2 , w — u + iv <E B. 
Of course ||W(u, u)|| =0 for u + iv ^ B. Therefore, we have 



(189) 



|W(u,«)|| 



~kn\v\ 3/2 

du dv < Kn I I — — ■ — ^— ; -r du dv 



|(u + iv) — (x + iy)\ J J B \(u + iv) — (x + iy)\ 

rd+25 , 

-kn\v\ 3 ' 2 [ ^ 



S r(3+28 

= Knl \v\e- kn ^ /2 = dv. 



-8 Ja-28 y/{u — x) 2 + (v — y) 2 

Now we will show that there is a constant C > such that 

f p+2S du ( 1 \ 

(190) I(x,y,v):= < Clog l + ] , , (x,y)eR 2 , veR. 

J a -28 y/(u- x) 2 + (v -y) 2 \ \v-y\J 

Indeed, since \{u + iv) — (x + iy)\ > \v — y\, on the one hand we have 

J a -28 \v-y\ \v-y\ \v-y\ 

We will use this estimate when \v — y\ is large. On the other hand, for \v — y\ small we have the following. 
Firstly, since \{u + iv) — (x + iy)\ > \u — x\, 

f p+2S du f^ 28 

(192) x<a- 2(5-1 or x>/3 + 25+1 =J> I{x,y,v) < f ] n</ du = f3-a + 4S. 

Secondly, we have may evaluate I{x,y,v) explicitly as 

, (B + 25-x\ , (a-25-x 

(193) I(x,y,v) = arcsmh [ — j j — ) — arcsmh 

from which it follows that 
(194) 

ff3 + 25 — x\ f(3 — a + 4(5+1 

a — 2<5— l<x<a — 2(5</3 + 2<5 =>• I(x, y, v) < arcsinh ( — j j — ) < arcsinh 

25 



a-25 \U — X\ Ja-28 



\v-y\ J V \ v -y\ 



v-y\ J \ \v-y 



a -25 <x < 13 + 25 



(195) 



and 
(196) 

a -25 < (3 + 25 < x < (3 + 25+1 



I(x,y,v) — arcsinh 
< 2 arcsinh 

I(x, y, v) < arcsinh 



(3 + 25-x 

\v - y\ 

(3-a + A5 
\v - y\ 

x — a + 25 
\v - y\ 



arcsinh 



x — a + 25 
\v - y\ 



< arcsinh 



(3 - a + 25 + 1 
\v - y\ 



The estimates (192 1 and ( 194 1— ( 196 1 may be combined for \v — y\ sufficiently small to give 

1 



(197) 



I{x, y, v) < C 2 log 



\v - y\ 



This estimate is useful for \v — y\ small. Taking (1971 with (1911 yields (1901. 
Using (1901 in (1891 and extending the integration from \v\ < 5 to v € M yields 



(198) 



\W{u,v)\\ 



\(u + iv) — (x + iy)\ 



du dv < Kn 



+oo 



-fe»M 3/2 io g ( i 



\v - y\ 



dv 



for some modified positive constant K > independent of (x, y) € M 2 and n. Now rescaling the integration 
variable by v — n~ 2 / 3 s gives 

(199) 



/ + OO / i 



dv = n- 1/3 



+ 00 



2 



ls|e- fc|s| log 1+ 

+ C30 



2/3 



ds 



\s\e " fc|s|3/2 log ( n 2 / 3 - 



+oo 



n 



2/3 



\s — n 2 /^y\ 



ds 



Wlogln) / \s\e- k ^ ' ds 



OC- 



-V3 / | s | e - fc l s l 3/2 log 1 



\s — n 2 / 3 y| 



ds. 



The first integral is clearly finite and independent of n, and since by Cauchy-Schwarz, 



(200) 



+ 00 



s\e- k ^ 3/2 logfl 



\s — n 2 / 3 y\ 



ds < 



s 2 e~ 2k ^ 3/2 ds 



1/2 



log 1 



ds 



1/2 



the second integral is bounded by a finite quantity independent of n. This proves that 



(201) 



sup 



T/ V^hM du dv < cr n -l/3 lQ (n) 

\{u + iv) - (x + iy)\ ~ 6V ' 



holds for some constant C > independent of n, for n sufficiently large. 

The Neumann series for the integral equation (182) corresponding to d Problem [I] therefore converges in 
L°°(IR 2 ) for sufficiently large n, and the estimate (185) follows immediately. □ 



7. LARGE-n ASYMPTOTICS FOR A(z) AND THE ORTHOGONAL POLYNOMIALS 

We will restrict our attention to fi + in "the bulk" (i.e. away from the endpoints a and 0) and also to 
the upper half-plane in the vicinity of the endpoint [3. Considerations for z near a are nearly identical, 
and we will omit them for the sake of brevity. Since the orthogonal polynomials being considered have real 
coefficients, the asymptotic behavior for z in the lower half-plane, C_, may always be obtained by complex 
conjugation. In this section and the next section we will use the notation 

(202) A„ := n- 1/3 log(n). 
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7.1. Asymptotics of the leading coefficients. The leading coefficients K n -i,n-i > and K n ,n > are 

obtained from A(z) using (60 1. In a neighborhood of z — oo, the matrix A(z) is related to E(x, y) via 



(203) 



A(z) 



i /2 E(x,y)i>(z)e n9( - z),J3 e 



nla 3 /2 



Using ( |166| , ( [175] ), and ([185]), we get 

A ia (z) = I e "^*» Uz) - ^z)- 1 + O ( v) ) 



A 21 ( Z ) = -i e «(^M) (tW -T(^)- 1 +0 (^)) 



(204) 



where the error terms are valid for z near z = oo. Now, from (|46[) we have 
(205) g(z)^log(z) + 0(z- 1 ), z -> oo, 



and from (164) and the condition that 7(2) — > 1 as z —> 00 we easily obtain 

1 



(206) 



7 (z) 



I - — (0 - a) + 0(z- 2 ), z^oo. 



Therefore, using (60 1, we obtain 
(207) = ^< 



'(1 + 0(A„)) and k= 



e""* (1 + 0(A„)) 



as n — > 00. These formulae may be combined with the asymptotic formulae for An(z) and A 22 {z) to be 
given below to obtain asymptotic formulae for the orthonormal polynomials p„_i(z) and p n (z). 

7.2. Asymptotics of the orthogonal polynomials in the bulk. Within the set f2+, the solution A(z) 
to Riemann-Hilbert Problem [l] is related to E(x,y) via 



(208) 



A(z) 



as follows from (61 1, (1391, and (1781. According to (59 1, the first column of A(z) contains the orthogonal 
polynomials: 



(209) 



1 2ni 

An (2) = Pn(z), and A 2 i{z) = p n -i{z). 

K"n,n ^"n— l.n — 1 



Via (208 1 these may be expressed as follows: 

(210) A n (z) = [E n (x,y) (t> n (z) + e~ ine ^ D 12 (z)) 

+ E 12 (x, y) (D 21 (z) + e- me ^D 22 (z)) 

(211) A 21 (z) = [E 21 (x,y) (Z) n (z) + e~ me ^D 12 (z)) 

+ E 22 (x, y) (b 21 (z) + e' me ^D 22 (z) 



a n{g{z)-t) 



Using the asymptotic estimate (1851 of E(x, y) — I, the relation (1751 between D(z) and D(z), and the 
explicit formula (166) for D(z) valid for z G f)+, straightforward manipulations yield 

(212) Au(z) = e n{9{z) - ie{x - v)/2] a{z) 



cos( -(n&(x,y)-<p(z)) ) (1 + 0(A„)) 



+ sm[-{ne(x,y)+<p(z)))0{A n ) 



(213) 



A 21 (z) 



,n{g(z)-Z-i&(x,y)/2) 



a(z) 



sin( -(nB{x,y)+<p(z)) ) (l + 0(A n )) 



+ cos - (n6(x, y) - <p(z)) O (A„) 
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where 
(214) 



a(z) 



(z - a)V*Q3 - z)i/4 



and <p(z) := arcsin 



2z - (a + /?) 

/3 — a 



are both functions analytic for z £ C \ (M. \ [a,P]). One apparent difficulty with these asymptotic formulae 
is that they involve an extension <d(x,y) of 9{x) that is completely arbitrary except that it must satisfy 
Properties 1-3 of Lemma [T| (While our proof of Lemma [I] was by construction, there was no assertion of 
uniqueness, and indeed there are many extensions Q(x,y) having the required properties.) On the other 
hand, it is also easy to see that the differences between various extensions Q(x,y) may be absorbed into the 
error terms. For example, if we fix x € (a, (3) and fix y > sufficiently small so as to be in the region f2 + , 
then by Property 2 of Lemma[T]we have Im(0(x, y)) < — ky for some k > (here we are using the fact that 
x is bounded away from a and f3), so (2121 and ( |213[ ) become 



(215) 

and 

(216) 



Au(z) = l -e n ^a{z)e-^' 2 (1 + O (A„)) , y » n" 1 , 



A 21 (z) = -i e ™(^M) a (z) e ^)/ 2 (1 + O (A n )) , y » n" 1 . 



On the other hand, if we suppose that x € (a, /3) is fixed and y — 0(n 1 ), then using Property 1 of Lemma [T] 
and the Mean Value Theorem we have 



(217) 



e(x,y) = 0(x) +e y (x,Ci)y, o<b<y 



and using (21 1 to eliminate <d y in terms of Q x and 90 this becomes 

(218) G(x,y) = 6(x)+iG x (x,£ 1 )y-2idG(x,Zx)y, 0<^<y. 

Since x is bounded away from a and (3, Property 1 of Lemma [T] guarantees further that Q xy is continuous, 
so another application of the Mean Value Theorem gives 



(219) 



= 9'{x)y + G xy {x,&)yti- 

Finally, using Property 2 of Lemma [T] to control 90 and using the assumption that y = C(n _1 ) we find 
(220) 6(ai, y) = 9(x) + iO'{x)y + 0(n~ 2 ). 



It follows that (|212| and (|213| become 
(221) 



and 
(222) 



A n (z) 



A 2 i{z) 



n(g( Z )-i9(x)/2+9'(x)y/2) 



a(x) 



-(n6(x) + inO\x)y - ^(a:)) + O (A„) 



-le 



n(g(z)-e-i8(x)/2+e'{x)y/2) 



a(x) 



(n9(x) + in9'(x)y + <p(x)) + O (A„) 



In particular, for y — 0+ we can write these in the form 

(223) A U (X) = e n{cV{,x)+l-<t>{w))/2 a ^ cog Q ( n0 ^ ^ x ^\ + q ^ An6 n 



{cV{x)+l-<t>{x))/2 



and 
(224) 



A 21 (x) = -ie n( - cV ^- l ~^ x »/ 2 a{x)sm Q (n9(x) - cp(x))J + O (A„e" 



{cV{x)-t-4>(x))/2 



2S 



A 11 (z) = V^e n ^ [F r \(z) (1 + (A„)) + Fl(z)0 (A„)] , 
A 2l {z) = 0Fe"toW-fl [f 2 (z) (1 + O (A„)) + f*(z)0 (An)] , 



7.3. Asymptotics of the orthogonal polynomials at the edge. Next, suppose that z £ S 1 ^ (see Fig- 
ure]^]). Then similar calculations in which the formulae ( 167[ ) and (1701 are used to find D(z) yield 

(225) 

(226) 

where 

Fl{z) := n V*w(z) [ e -"/ 3 e n ^« 3/2 / 2 Ai ((f ) 2/3 U/5 (z)e 2 "/ 3 

(227) + e W3 e -nM*) 3/ 72 e -*ne(^) Ai ((3«) 2 / 3 ^e"^ 

+ n- 1 ^w(z)- 1 [ e -2-/3 e n^( Z ) 3/2 /2 Ai ' ((3«)2/3^ (z)e2W3 ^ 

+ e 2«/3 e -n^( 2 ) 3/2 /2 e -me(., a ) Ai / ((3«) 2 / 3 U/J ( z ) g" 2 " 4 / 3 ) 



F 2 (z) := 



-5ni/6 nu l3 (z) 3/2 /2 



Ai 



3n\ 2 /3 



U(3(z)e 



2ni/3 



(228) 



+ e -«/6 e -^( z ) 3/2 /2 e -me(x,s) Ai ((f ) 2/3 u p (z)e- 2 ^ 3 



and 
(229) 



w(z) 



5) 



1/6 



Now it will also be useful to have the asymptotic behavior of A n (z) and A 2 i(z) for z S ff| P I C + (see 
Figure [4|, and for this purpose we note that for such z we have B(z) = D(z), so in place of (2081 we have 
instead 

(230) A(z) = e nea3/2 E(x, y)t>{z)e n9 ^ e- nta3 / 2 , z = z + iy E C \ (fi + U U ft a U Q ). 
Therefore, for such z: 

(231) A u (z) = \E n {x,y)D n {z) + E 12 {x,y)D 21 (z)\ e n ^ z \ 

(232) A 21 (z) = \E 21 (x,y)D 11 (z) + E 22 (x,y)D 21 (z<\ *»C»(*M) 



Supposing that z £ Sjfl C+, we may now proceed by using ( [T67| and ( [169} to find D(z), with the result 
that 



A n (z) = 0Fe B »Me™*« 8/a / 2 [ G i ( Z ) (1 + ( A „)) + Gl(z)0 (A„)] , 
A 21 (z) = ^ e «(^)^) e »M*) 3/2 /2 [G 2 (z) (1 + O (An)) + (z)O (A„)] , 

G,\(z) := n^^Ai ((f ) 2/3 U/3 (z)) - n-V^^J-iAi' ((f ) 2/3 U/3 (*)) , 
G 2 (z) :=-zn 1 /6 w(z )Ai(( 3 f) 2/3 M/3 (z)) - m" 1 / 6 ^)" 1 ^' ((f ) 2/ %(z) 



(233) 
(234) 
where 
(235) 

(236) 

If we assume that C : — (3n/4) 2 / 3 up(z) is bounded, then Property 3 of Lemma [I] guarantees that 
e - m e(x : y) = e 2C 3/2 /3( 1 + 0( n -2/3^ and so with the use of (TfygJ we see that ( gggj ) an d p26| agree, re- 
spectively, with (2331 and (2341 up to error terms; we therefore have 

^3/2 , 



(237) 



Au(z) 



fiHwtfi)nV*(X /3 Ai(C) + O [n^ e A n 
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and 
(238) 



A 21 (z) 



n(g(z)-() 



-i^(/3)n 1/6 e 2 ^ 3/2 / 3 Ai(C) + O (n 1/6 A„)~ 

for £ bounded with < arg(£) < tt, where 
(239) 

Moreover, we may observe that for z near (3 , and z G C+, the following local expansion holds true: 
(240) 



= /3 + (An)- 2 / 3 C, A:= i [-^G9)]" 1 . 



(This is proved in the Appendix under the conditions on the external field V in force in this paper.) Therefore, 
(237JI and (2381 may also be written as 



(241) 

and 
(242) 



A n (/5 + (An)- 2 / 3 c) 



A/6„n(cV(f3)+£)/2 fl n 1 ' 3 cV'(l3)\- 2 ' 3 <;/2^ 



O 



(nV«A n « 



n( C y(/3)+<')/2 p n 1 / 3 cy'(/3)A- 2 / 3 C/2 



A 21 (/3 + (An)- 2 / 3 c) = -m 1 / 6 e n ( cy W-^/ 2 e" 1/3cV '(« A " :!/3 «/ 2 V^i ( ;( / 9)Ai(C) 

+ O ^ n 1 /6 Arie «(cV(/3)-£)/2 e n 1 /3 c ^ (/ 3 )A -2/3 c/2 X 



Here we used the fact that the error term in (240) is 0(n 4 / 3 ) for \(\ bounded, so the dominant terms in 
the errors come from ( 237 1 and ( 238 1 . 



8. ASYMPTOTICS FOR THE DERIVATIVE AND APPLICATIONS TO RANDOM MATRIX THEORY 

We have established large-n asymptotics uniform with respect to z for the matrix A(z), for all z S C. In 
particular, An(z) and A 2 i(z) possess asymptotic descriptions in a neighborhood of the interval [a,0\. In 
this section we derive asymptotic descriptions for the derivatives A' n (z) and A' 2 i(z), both in "the bulk", i.e. 
for x £ {a, 13) as well as near the endpoints a and f3. (In fact we will only consider the endpoint /?.) 

Our aim is to obtain derivative asymptotics in order to establish bulk and edge universality for unitar- 
ily invariant matrix models with external fields that possess only two Lipschitz continuous derivatives, as 
described in the Introduction. It is by now well-known (see, for example, |16j ) that if one obtains an as- 
ymptotic description of the orthogonal polynomials and their derivatives of the form which we obtain here, 
then the corresponding asymptotic formulae for the reproducing kernels follows and exhibits universality 
(independence of details of the external field V), and so we will omit these details. 

8.1. Analysis of derivatives in the bulk. Let x € (a, f3) be bounded away from the endpoints as n — * oo. 
Since A\\(z) and A 2 i(z) are polynomials and hence entire functions, we may express their derivatives at x 
by Cauchy's integral formula: 



(243) 



A' n {x) 



1 

2~7Ti 



Aji(s) ds 
(s — x) 2 ' 



J = 1,2. 



Moreover, since (see ( 209 1 ) An(z) and iA 2 i(z) have real coefficients, we may use complex-conjugation 

A 2 i(s) ds* 



symmetry and the reality of x to write 

77 \Jr \ s ~ x ) 



(244) 



and A 2 i(x) = — :Re 



where T is any path of integration that begins on the real axis to the right of x and terminates on the real 
axis to the left of x, and that avoids the singularity at x by passing through the upper half- plane. For our 
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calculations, we will take the path T to be a semicircle of radius n 1 centered at s — x: s = x + n 1 e luJ for 
< to < 7r. Thus we have 



A' n (x) = ^Im (J A n (x + n- l e iul )ie- iui doj\ , 
A' 21 {x) = ^Rc (J A 21 (x + n- 1 e luJ )ie-^ duj 



(245) 
(246) 

We will now substitute from (2121 and (213), but first we write them in a more suitable form. Since a(z) 
is bounded for z £ T, and since g(x + iy), &(x,y), and ip(x + iy) are all differentiable, Taylor expansion 
about x — shows that for z E T we have both 

(247) A u (z) = hh{x, y) + (A„e" ff +^) and e ni A 21 (z) = M 2 (x, y) + O (& n e n ^) , 
where g+{x) is the boundary value taken by g(z) as z — > x from C+, and where {z = x + iy) 

(248) Mx(x, y) := e n{a{z) " lB{x ' v)/2) a{z) cos Q(n6(x, y) ~ ip(z))j 



and 
(249) 



(x,y) := -ze n(9(2) - ie(a; ^ /2) a(z)sin (^(nQ(x, y) + <p(z)) \ . 



One important observation is that Mi(x, y) and iM 2 (x, y) have real boundary values taken on the real axis 
from the upper half-plane. Indeed, the analytic functions a(z) and f(z) are real for real z, Property 1 of 
Lemma [T] implies Q(x, 0) = 6{x) S JR., and furthermore by ( |49| , 

(250) 2g + (x) - i6(x, 0) = 2g + (x) — i6(x) = g+{x) + g~(x) G E, a: G (a,/3). 



Using ( |247j in ( [245J and ( [246] ) gives 

(251) A' n (x) = ^Im (J Mi(x + n' 1 cos(w), sin(w))ie" 1 " + O (nA„e" ff + (:E) ) , 

(252) e^jfc) = ^.Re Qf M 2 (:z + n' 1 cos(w), n" 1 sia(w))ie- <w + O (nA n e n9+{x) ^ . 

We begin our analysis by integrating by parts: since for M = Mi or M = M 2 , 

r d 
M(x + n' 1 cos(u),n~ l siYL(<jj))ie~ lu) duj = - / M(x + n^ 1 cos(w), n~ x sin(w)) — (e~ luJ ) 



(253) 



Jo 

M(sc - n _1 , 0) + M(x + n~ l , 0) 
1 
/? 
1 



(i 



H — / M y (x + n 1 cos(w),n 1 sin(w)) cos(w)e luJ du> 
M x (x + n^ 1 cos(cj), nT 1 sin(w)) sin(a;)e~ l " duj, 



o 



we see that upon taking the imaginary part (for A'^x)) or real part (for A' 21 (x)), the boundary terms vanish: 
(254) nlm ( [ Mi(x + n" 1 cos(w), sin(w))ie~ lw du^ 

Im ( / Mi y (x + n" 1 cos(oj), n _1 sin(w)) cos(w)e~ lw cLj 



Im 



Mi x (x + n 1 cos(cj),n 1 sin(w)) sin(o;)e lu duj I , 
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and 



(255) nRefj M 2 (x + rT 1 cos(u), n~ x sin(w))ie" iw dw 

M2 V (x + n cos(w), n~ l sin(w)) cos(w)e _l " dui 



Re 



Re 



u - - \ 

M 2x {x + n~ x cos(w), n~ x sin(w)) sin(a;)e _l " dw ) . 



Moreover, using (211 to eliminate Mj y in favor of Mj X and dMj for j = 1,2, these become 



(256) nlm 



Mi(;c + n 1 cos(w),n 1 sin(u;))ie lw e?w 
= Rc 



Mi^rr + n 1 cos(a->),n 1 sin(w)) g?oj 
2Re ( / dM x (x + n~ l cos(w), n" 1 sin(w)) cos(uj)e~' lu] duj I , 



M 2 (x + n cos(u;),7i 1 sin(w))ie lw duj 
= -Im 



Ma^x + n 1 cos(cj),?i 1 sin(cj)) e&J 

dM2(x + n" 1 cos(w), n" 1 sin(u)) cos(w)e _ia ' du 



and 

(257) nRc 



+ 2Im 

By splitting cosines and sines into exponentials, we may write M\{x, y) and M2{x,y) in the form 

(258) Mi(aj,y) := u-{z)e na{z) + u + {z)e n{g{z) - m{x > v)) 
and 

(259) M 2 (x,y) := u~ {z)e n ^- i@( ' x ' y ^ - u+{z)e n9(z \ 
where 

(260) 



u ± (2) := ia(z)e ±l¥, W/ 2 . 



One reason for writing Mi(x,y) and M2{x,y) in this way is to explicitly display their dependence on n; 
indeed u ± (z), g(z), and 0(a;,y) are independent of n. Using the fact that u ± (z) and g(z) are all analytic 
functions of z — x + iy on the contour T, we have 
(261) 

flMi(a?,y) = -mu+(z)e"( 9 ( z ^ te ( a; ^»9e(x,y) and dM 2 (x,y) = -mu-(z)e n ( ! '( 2 )- ,e < I ' ! '»»9e(i,i/). 

Now, u ± (z) are bounded functions, and according to Property 2 of Lemma [I] d<d(x,y) = 0(n _1 ) for 
x + iy G r. Since also g(z) and Q(x,y) are differentiable and (by Property 1 of Lemma[T]) Q(x,0) — 9{x) G M, 
we finally learn that 

(262) dMx(x, y) = (e" 9+w ) and dM 2 (x, y) = O (e^+^A , x + iyeT. 
Therefore, (251) and (252) may be written in the form 

(263) A' n (x) = ^Rc Qf M lx (x + n- 1 cosiuj),^ 1 sm(oj)) du?j +0 (e ng +^ 
and 

(264) e n£ A' 21 (x) = ~Im (J M 2x (x + rC 1 cos(w), rT 1 sin(o;)) du?j + O (e ng +^ . 
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Now, each term in Mi(x,y) and M2(x,y) is of the form u ± {z)e n ^ 9 ^ la ®( x >y)) where either a — or a = 1, 
and 

9 



(265) 



Ox 



u (z)e 



n(g{z)-ia&{x,y)) 



(z) + rm±(z) (g'(z) - *a6 a; ( 2 :, y))] e "0fW-^(*.v)), 



so evaluating for z = a; + n~ 1 e M we may expand the result for large n. Using Property 1 of Lemma [T] to 
assert second-order differentiability of Q(x,y) in the bulk and Property 2 of the same Lemma to guarantee 
that d<d(x,0) — 0, we finally arrive at 
(266) 

^ [ u ±( 2 ) e ™(9«-^e(=^))' = nu ±( x ) (g' + (x) - ia6'(x)) e <9+{x)-i«e{x)) ^e™ + e> (e" 9 



ig+{x) 



z—x+n^ 1 e ruJ 

where q := g' + (x) — ia9'(x) is independent of lj. Since for any complex number q 

(267) f e qei " duj = n, 

Jo 

we obtain from ( p63| and ( p64| that 

A' u (a) = Re (n U - (x)g' + {x)e n9 + {x) + nu+(x) (g' + (x) - id'(x)) e »(ff+(*)-*?(*)A + Q L,ng+{x)\ 



(268) 



Re 

d 

dx 



d 

dx 



M 1 (x,0+) \ +o(e ng +( x A 



McV(x)+e-<p(x))/2 a ^ C0S n {n6{x) _ v ^ 



O 



(f 



(cV(x)+l-<P(x))/2 



and 



A' 21 (x) = flm (nu~(x) (g' + (x) - iB'{x)) e M9+(x)-i6(x)-t) _ nu +( x ) g i + ( x ) e n(g+(x)-^ + ^ e n{g+(x)-lj^ 



(269) 



= ilm I ( 

d 

dx 



— M 2 (x, 0+) j + O (e" (9 + (a:) -^) 

le n(cV(x)-i-^))/2 a{x) sin Q(„0( x ) + ^(^ + O ^(cV^M-fl*))/^ 



Comparing these results with ( 223 ) and ( 224 1 shows that the asymptotic formulae for the derivatives of the 
orthogonal polynomials on the real axis in the bulk may be obtained from the corresponding asymptotic 
formulae for the polynomials themselves by differentiating the leading terms. 

8.2. Analysis of derivatives at the edge. We will now apply similar considerations to the asymptotic 
formulae ( 237 1— ( 238 1 , to obtain asymptotics for derivatives of the orthogonal polynomials that are valid for 
z in a vicinity of the endpoints (3 and a. We will present the details for the cndpoint z = (3, as the argument 
for the behavior near z — a is entirely similar. More precisely, our aim is to establish asymptotic formulae 
for the quantities 



(270) 



d 

dC 



A n (/3+(An)- 2 / 3 c) 



and 



d( 



A 2 i(/3+(An)~ 2 / 3 c), A 



which arc what one needs to establish universality of the distribution of the largest eigenvalue in Hcrmitian 
random matrix theory. 
Let 



(271) 
(272) 



n ( C ) ;= e -n(cV { p )+ l)/2 e -n^icV>{P)X-i'\/2 Aii U + {Xn)~^\) 



MO 



._ -n(cV(0)- 



According to (241 1 and (2421, these may be expressed as 

(273) n (C) = n^V^u^AitC) + O fn 1/6 A„ 



:!3 



and 
(274) 



75,(0 = -m 1/6 ^(/3)Ai(C) + O (n x / 6 A n ) 



for \(\ bounded. 

Since Ti(£) and T2(C) axe entire functions that are real for real £> just as in the analysis in the bulk we 
may express the derivatives t[(Q and t' 2 (£) in terms of Cauchy's formula as 



t{(0 = -Im 

7T 



7-^(0 = ^Re ( [ 



r 2 (0^ 



(275) 
and 

as long as £ € M, where T is any path in the upper half-plane from the real axis to the right of £ to another 
point on the real axis to the left of (. But, since the dominant terms in (2731 and (2741 are entire functions 
of C, and since |£ — £| is bounded away from zero on the contour T of finite length, it follows from a residue 
calculation that 

d 



(277) 
and 
(278) 

Therefore, 



r{(0 = 



4(0 



d_ 



(-zn 1 / 6 V^«;(/3)Ai(C)) + O (n 1 / 6 A„) 



(279) 



(280) 



d^ 



d< 



-A 



ii (/3+(An)- 2 / 3 c) = ^ 1 / 2 c r( / 3)A- 2 / 3 e"( cV '(«^/ 2 e nl/3cy '(' 3 ) A " /3 «/ 2 ^(/3)Ai(C) 
+ O rn 1 / 2 A n e" (cVr(/3)+ ^ /2 p nl/3cV " (/3)A_3/3c/2 " 

A 21 (/3+(An)- 2 / 3 c) 



^^n 1 / 2 cT/'(/3)A- 2 / 3 e"( cV ^^/ 2 e " 1/3cl/ '^ A " /3<; / 2 ^(/3)Ai(0 



These may also be written as 
d 



(281) 



(282) 



dC 
d 



O (n 1 / 2 A n e n 

d 



(cy(/3)-£)/2 e n 1 / 3 C y'(/3)A- 2 / 3 C/2 



An (p+ (An)- 2 / 3 C) = £ »V6 e »<cV(fl^)/V 1/ScV W- a/ Va^( /3 ) Ai (0 



A 21 (/3+(An)- 2 / 3 C 



-1 

d 



(n 1/2 A„e™ (cy(/3)+£)/2 e" 1/3cy ' (/3)A_2/3c/2 ) 
-zn 1 / 6 e"( cV («- £ )/ 2 e" 1/3cV "(« A " /3 ^ 2 ^(/3)Ai(C) 



+ C (n 1/2 A„e" (cy(/3) ^ )/2 e™ 1/3cV " (/3)A_2/3c/2 ^) , 



which, upon comparing with (241 1 and (242 1, show that asymptotic formulae for derivatives valid at the 
edge may be obtained by differentiating the leading terms of the corresponding formulae for the polynomials 
themselves. 

Appendix: Convex External Fields 

Everywhere in this Appendix we shall assume that (i) the external field grows sufficiently rapidly as 
\x\ — > co, and (ii) the external field is strictly convex, and possesses d continuous derivatives, with d > 2. 
In a separate discussion below, we will consider the specific situation that V possesses just two Lipschitz 
continuous derivatives. The main results are summarized in Lemma [3] at the end of this Appendix. 

The assumed growth and strict convexity of V(x) and the positivity of c implies that /it* is compactly 
supported and absolutely continuous with respect to Lebesgue measure, with support consisting of a single 
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interval [a, (3\ for some real a < (3. To obtain a formula for /i* in this case, we consider the auxiliary function 
g(z) defined in ( ]46[ ). analytic for z G C \ (—00, /3] (in the present case the integral is taken over the single 
interval [a,/3]). In terms of g(z) the variational condition (44) becomes ( 47 1 which we rewrite here: 



(283) 



cV(x) - (g+(x) + g-(x)) — —I, a < x < (3, 



where g+{x) and g~(x) denote the boundary values taken by g{z) as z — » x with z G C + and z G 
respectively. 

Assuming that differentiation commutes with taking boundary values (47 1 and (48 1 imply that 

g' + {x) + g'_{x) 



(284) 



g'+{x) -gL(x) 



cV'(x), 
0, 



x G (a, (3) 
x 6l\ (a, 0). 



In particular, g'(z) is an analytic function for z G C \ [a,/3]. To find g'(z) from these conditions, we 
introduce the function R(z) satisfying R(z) 2 = (z — a)(z — (3) such that R(z) is analytic for z G C \ \ot,f3] 
and R{z) = z + 0(1) as z -> 00. Setting 3' (2) = f(z)R(z) for some new unknown function f(z), we find 
that like g'(z) and R(z), /(z) is an analytic function of z for z G C \ [a,/?], and that its boundary values 
taken on (a, (3) from the upper and lower half-planes satisfy the relation 



(285) 



/+ Or) -/_(*) 



cV'(x) 

R+{x)' 



a < x < (3 . 



Since g'(z) = l/z + 0(l/z 2 ) as z — > 00, it follows that f(z) — 1/z 2 + 0(l/z 3 ) as z — > 00, and hence 



(286) 



/(*) 



V 7 (a) 

2^y Q (s-z)R+(sy 



Considering (286 1 for large z, we sec that 

* V'(s)ds 



(287) 



R + (s) 



= 0, 



^ sV'(s) ds 



= -2m. 



These two equations determine the endpoints a and f3. With a and (3 chosen so that the equations ( |287| 
hold, we may obtain a formula, in terms of a Cauchy principal value integral, for the density ip( x ) of the 
equilibrium measure valid in the support interval a < x < (3: 

cR+{x) [ f) V'(s)ds 



(288) 



J a. 



2t 2 J a (s-x)R+(s)' 
Defining a real-valued function h(x) for x G K by the formula 



(289) h(x) := - 

7T 

it is straightforward to verify that 
(290) 
and that 



V'(s) - V'(x) ds 



tp(x) = —R+(x)h(x), 
2m 



(291) 



(a?) := cV(ar) + 1 - g+(x) - g-{x) = < 



R + ( s y 

a < x < (3, 

R(s)h(s) ds, 
R(s)h(s) ds, 





x < a 



x > (3. 



From the assumption that V(x) is d times continuously differentiable, we see that h{x) is d — 2 times 
continuously differentiable. Also, since i/R + (x) is positive for a < x < [3 and 



(292) 



V» - V'(z) 



V"(ts + (l-t)z)dt, 



the assumption of convexity of V(x) implies that h(x) is strictly positive for all Also, since g+(/3) — 

<7_(/3) = and since for a < x < (3 we have g' + {x) — g'-(x) = —2iriip(x), we conclude that 

r fi 

(293) 9{x) := -i(g+(x) - g-{x)) = -ic / R+(s)h(s) ds, a < x < (3. 



We may now use (291 1 to verify (from the positivity of h{x) and the facts that R{x) > for x > (3 while 
R(x) < for x < a) the strict inequality 4>{x) > for x < a and x > (3, as required by Condition [2] Also, 
from ( |293[ ) we see that < 9(x) < 2tt and tp(x) > (i.e. 9'(x) < 0) both hold strictly for a < x < (3 as 
required by Condition [2j Moreover, since h(s) is d — 2 times continuously differentiable, 9{x) is d — 1 times 
continuously differentiable for a < x < [3. 

One may prove ( |240[ ) as follows. From the identity g'{z) = R(z)f(z), we have that for z near /?, 



(294) /(z) = 



cV'((3) 
2R(z) 



c 
2~Tri 



0V'{a)-V'((3) dfj l cV"(/3) 



(s-/3)R + (s) 



2R(z) 



Zm 



-{z-(3) 



f V'(s) - V'(f3) 

8-/3 



V"((3) 



ds 



(s-z)R + (s) 



Recalling that 2g((3) — cV(/3) — I — 0, we learn that 



(295) g(z) 



cV({3) 



cV'(P) 



(z-[3) + 




V'(s) - V>((3) 



(s-f3)R+(s) 
c{z' - f3)R{z') r 13 (V'{a)-V'{p) 




R(s)ds+ cV ^(z-(3) 2 



2ni 







V"{p) 



ds 

(s-z')R+( S ) 



dz'. 



Now, recalling that 

(296) 
and 
(297) 



^) 3/2 —^(-/3) 3/2 , 



hp(x) 



we may use (295) to obtain the value of h'p((3): 

2c4W~- 



d i x ) - 9+{x)-g-{x) 
— — i t=^= , a < x < p, 



(298) 



h'J(3) 
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13 V'(s) ~ V'((3) 
: (s-P)R + (s) 



ds < 0, 



and then the expansion ( 240 1 follows by adding ( 295 1 and ( 296 ) . 



Regularity for V" Lipschitz. It is also straightforward to derive a formula for g"(z). This is a useful 
exercise if one assumes only that V" is Lipschitz continuous, which we do throughout this subsection. 



By further differentiation, (47 1 and (48 1 imply that 

cV"(x), 



(299) 



g'i(x)+g'L{x) 



x £ (a, (3) 
x E 



g'{{x) - g"_{x) = 0, i6l\(a,j3), 

In particular, g"{z) is an analytic function for z € C \ [a, 0\. To find g"(z) from these conditions, we set 
g"(z) — F(z)/R(z) for some new unknown function F(z), we find that like g"(z) and R(z), F(z) is an 
analytic function of z for z € C \ and that its boundary values taken on (a,/3) from the upper and 

lower half-planes satisfy the relation 

(300) F+(x) -F_(x) =cV"(x)R+(x), a < x < (3 . 
It follows that F(z) is a function of the form 

(301) F(z) 



c 
2t?l 



V ,, (s)R + (s)ds 
(s-z) 



:;(> 



Since ijj'(x) = — 2ni(g'l (x) — g'!_{x)) for a < x < (3, we obtain 

2wi 



ip'(x) 



(302) 



R+(x) 
2c 



(F + (x)+F_(x)) 



P V"(s) - V"{x) 2nicV"(x) 

R + {s) ds + 

s — x R+{x) 



a + (3 



from which it follows that y/(x — a) {(3 — x)ip' {x) is bounded uniformly for a < x < /3. In other words, 

C 



(303) 



\ip'(x)\ < a<x</3. 
V(z - a)(f3 - x) 



From the above considerations we have the following formula, which is valid for x < a and also for x > (3: 



(304) 



<P"(x) 



ic 



ttR(x) J a 



V"( S )-V"(x) cV"(x) 

R + {s)ds + 

s — x R(x) 



a + f3 



which in turn implies that \x — a] 1 ^ 2 ^ — I3\ 1 / 2 (j)" {x) is bounded on any compact subset of K. (Of course, 
this quantity diverges as \x\ — > oo, with x e M.) 

On the other hand, one may also consider the quantity F(z) := g{z)/ R(z), for which the following identity 
can be shown to hold true: 

ds \ c r p V{s)-V{z) ds ^cV(z) + £ 
,R(s)(s-z)J 2iriJ a s-z 



(305) 



F±(z) 



One direct consequence of this last identity is that 

ds 



(306) 



9+{z) -9-{z) = R+(z) 



-oo R(s)(s - z) ni 



R+(s) 2R ± {z) 

13 V(s) - V(z) ds 
, s-z R + (s) 



z e 



, ze(a,f3], 



with both quantities appearing within the square brackets on the right hand side of ( 306 1 possessing at least 
one Lipschitz continuous derivative for all z € [a + e, f3] for any e > 0. Similarly, we have 



(307) g+(z) + g.(z)-cV(z)-l = R(z) 



ds 



R(s)(s — z) 7ri J a 



V(s) - V(z) ds 
s — z R+(s) 



z € [/3, oo), 



and again the quantities appearing within the square brackets on the right hand side of (307) possess at 
least one Lipschitz continuous derivative for all z S [f3, oo). 

The behavior near z — a is slightly more subtle, but using the identity 



(308) 



ds 



R{s)(s~z) R ( z ) J^Ris^s-z) 



ds 



where R(z) — sgn (Im(z)) R(z) is the function which coincides with R(z) in C+, and is analytic in C \ 
((— oo,a] U [f3, oo)). Indeed, using (3081, the identity (3061 becomes 



(309) 



j+(z)-g-(z) = 2m + R+(z) 



ds 



V(s) - V(z) ds 
R(s)(s — z) iri J a s — z R+(s) 



, z £ [a, (3), 



and once again the quantity within the square brackets on the right-hand side of |309| possesses at least one 
Lipschitz continuous derivative. Similarly, the identity (3071 can be rewritten, in light of 308 as follows: 



(310) g+(z) + g_(z) - cV(z) - I = R(z) 



p ds c r p V(s) - V(z) ds 

R(s)(s — z) ni J a s — z R+(s) 



z 6 (— oo, a] 



the quantity within the square brackets again possessing one Lipschitz continuous derivative. 
We summarize the results of this Appendix with the following Lemma. 

Lemma 3. Suppose that the external field V possesses two Lipschitz continuous derivatives, is strictly convex, 
and grows faster than [log(l + x 2 )] 1+e for some e > 0. Then the density tp( x ) °f the equilibrium measure /i* 
is supported on a single interval, [ot,/3]. On this interval, the function ip has the following properties. 
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The function ip may be expressed via (2901 with h, defined in (2891, being Lipschitz continuous on 
[a, ft. 



• The function tp has one derivative, which satisfies the bound (303) 



In vicinities of the endpoints f3 and a, the related function 8 (recall 9'{x) — —2mp(x)) satisfies (cf. 



(3061 and (309)j 



(311) e( X )=2.[% {s)d s=h R ^y?> ^(a + e,/?) 

Jx [27r + iR + (x)h a (x), x G [ol,P — e), 

for some small e > 0, with hp and h a being positive functions on (a, (3). In addition, hp possesses 
one Lipschitz continuous derivative on (a,/3], and h a possesses one Lipschitz continuous derivative 
on [a,P). 

On the complementary set (-co, a) U oo), the following properties hold true: 

• The quantity 4>{x) :— cV{x) + I — g+(x) — 9-{x) possess two derivatives. The first derivative <fi'(x) 
may be obtained from ( 291| , with h(x) being Lipschitz continuous on (—00,0] U [/3,oo). 



The second derivative <^"(a;) satisfies the inequality 

(312) \x-a\ 1/2 \x- I3\ 1/2 \(f>"(x)\ < C 

on any compact subset ofM.. 
• In vicinities of the endpoints a and j3, the function 4>{x) satisfies (cf. (307) and (310 \) 

—R(x)hp(x), x€[/3, 00) 
R(x)h a (x), a; € (— 00, a], 



(313) <p(x) 



with hp and h a being extensions, to [j3, 00) and (—00, a] respectively, of the functions of the same 
name, formerly defined on and [a, (3), respectively. These extensions possess one Lipschitz 

continuous derivative as well. The function hp is strictly negative on all o/(/3, 00), and the function 
h a is strictly positive on all of (—00, a). 

Note that the functions h a {x) and hp{x) used in the main text are simply related to h a (x) and hp(x) as 
follows: 



(314) h a (x) — \J f3 — x h a {x) and hp(x) = \/x — ahp(x). 



p: 

[9 

[10: 
[11 
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